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Abstract 

This work continues our studies of nonlinear evolution of a system of wavepackets. 
We study a wave propagation governed by a nonlinear system of hyperbolic PDE's 
with constant coefficients with the initial data being a multi-wavepacket. By definition 
a general wavepacket has a well defined principal wave vector, and, as we proved in 
previous works, the nonlinear dynamics preserves systems of wavepackets and their 
principal wave vectors. Here we study the nonlinear evolution of a special class of 
wavepackets, namely particle-like wavepackets. A particle-like wavepacket is of a dual 
nature: on one hand, it is a wave with a well defined principal wave vector, on the 
other hand, it is a particle in the sense that it can be assigned a well defined position 
in the space. We prove that under the nonlinear evolution a generic multi-particle 
wavepacket remains to be a multi-particle wavepacket with a high accuracy, and every 
constituting single particle-like wavepacket not only preserves its principal wave number 
but also it has a well-defined space position evolving with a constant velocity which 
is its group velocity. Remarkably the described properties hold though the involved 
single particle-like wavepackets undergo nonlinear interactions and multiple collisions 
in the space. We also prove that if principal wavevectors of multi-particle wavepacket 
are generic, the result of nonlinear interactions between different wavepackets is small 
and the approximate linear superposition principle holds uniformly with respect to the 
initial spatial positions of wavepackets. 

1 Introduction 

The principal object of our studies here is a general nonlinear evolutionary system which 
describes wave propagation in homogeneous media governed by a hyperbolic PDE's in W^, 
d = 1,2,3,... is the space dimension, of the form 

5.U = --L(-iV)U + F(U), U(r,r)|^^o = h(r), r e < (1.1) 

Q 
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where (i) U = U (r, r), r G M*^, U G C^-^ is a 2J dimensional vector; (ii) L (— iV) is a 
linear self-adjoint differential (pseudodifferential) operator with constant coefficients with 
the symbol L(k), which is a Hermitian 2 J x 2 J matrix; (iii) F is a general polynomial 
nonlinearity; (iv) > is a small parameter. The properties of the linear part are described 
in terms of dispersion relations (k) (eigenvalues of the matrix L(k)). The form of the 
equation suggests that the processes described by it involve two time scales. Since the 
nonlinearity F (U) is of order one, nonlinear effects occur at times r of order one, whereas 
the natural time scale of linear effects, governed by the operator L with the coefficient 1/g, 
is of order g. Consequently, the small parameter g measures the ratio of the slow (nonlinear 
effects) time scale and the fast (linear effects) time scale. A typical example of an equation 
of the form fll.ip is the nonlinear Schrodinger equation (NLS) or a system of NLS's. Many 
more examples including a general nonlinear wave equation and the Maxwell equations in 
periodic media truncated to a finite number of bands are considered in [8], [9]. 

As in our previous works [8], [9] we consider here the nonlinear evolutionary system 
(11. ip with the initial data h (r) being a sum of wavepackets. The special focus of this paper 
is particle-like localized wavepackets which can be viewed as quasi-particles. Recall that a 
general wavepacket is defined as such a function h (r) that its Fourier transform h (k) is 
localized in /9-neighborhood of a single wavevector k^, called principal wavevector, where (3 
is a small parameter. The simplest example of a wavepacket is a function of the form 

h(/3;k)=/3-V'''-/i(^^^^g„(k,), kGM^ (1.2) 

where g„ (k^,) is an eigenvector of the matrix L (k*) and h (k) is a scalar Schwartz function 
(i.e. it is infinitely smooth and rapidly decaying one). Note that for h {f3, k) of the form (11. 2p 
we have its inverse Fourier transform 

h (/3; r) = h{P{r- r,)) e^''* )g„ (k,) , r G R'^. (1.3) 

Evidently, h (/3, r) described by the above formula is a plane wave e^''*'"g„ (k^.) modulated 
by a slowly varying amplitude (/3 (r — r*)) obtained from h (z) by a spatial shift along the 
vector r* with a subsequent dilation with a large factor ^. Clearly, the resulting amplitude 

has a typical spatial extension proportional to and the spatial shift produces a noticeable 
effect if |r*| ^ f3~^. The spatial form of the wavepacket (11.31) naturally allows to interpret 
G M'^ as its position and, consequently, to consider the wavepacket as a particle-like one 
with the position r,, G R'^. But how one can define a position for a general wavepacket? 
Note that that not every wavepacket is a particle-like one. For example, let, as before, the 
function h (r) be a scalar Schwartz function and let us consider a slightly more general than 
(II. 3p function 

h (/?; r) = [h {(3 (r - r,i)) + /i (/? (r - r.a))] e**^g„ (k,) , r G R^ (1.4) 

where r^,i and r*2 are two arbitrary, independent vector variables. The wave h (/?, r) defined 
by (ll.4p is a wavepacket with the wave number k* for any choice of vectors r.,,! and r=i,2, but 
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it is not a particle-like wavepacket, since it does not have a single position r*, but rather it 
is a sum of two particle-like wavepackets with two positions r^,i and r^,2- 

Our way to introduce a general particle-like wavepacket h k,,,, r*o; r) with a position 
r=KO is by treating it as a single element of a family of wavepackets h k*, r,,; r) with r* G M'^ 
being another independent parameter. In fact, we define the entire family of wavepackets 
h k^,, r^,; r), G M"' subject to certain conditions allowing to interpret any fixed G M'^ 
as the position of h k^., r^,; r). Since we would like, of course, a wavepacket to maintain 
its particle-like property under the nonlinear evolution, it is clear that its definition must 
be sufficiently flexible to accommodate the wavepacket evolutionary variations. In light 
of the above discussion the deflnition of the particle-like wavepacket with a transparent 
interpretation of its particle properties turns into the key element of the entire construction. 
It turns out that there is a precise description of a particle-like wavepacket, which is rather 
simple and physically transparent and such a description is provided in Deflnition 12.21 below, 
see also Remarks 12. 4[ 12.51 The concept of the position is applicable to very general functions, 
it does not require a parametrization of the whole family of solutions, which was used, for 
example in [32], |25], [26] . 

As in our previous works we are interested in nonlinear evolution not only a single particle- 
like wavepacket h k*, r^.; r) but a system {h k^;, r^,;; r)} of particle-like wavepackets 
which we call multi-particle wavepacket. Under certain natural conditions of genericity on k*^ 
we prove here that under the nonlinear evolution: (i) the multi-particle wavepacket remains 
to be a multi-particle wavepacket; (ii) the principal wavevectors k,.; remain constant; (ii) the 
spatial position r,,; of the corresponding wavepacket evolves with the constant velocity which 
is exactly its group velocity ^VcUn (k*i). The evolution of positions of wavepackets becomes 
the most simple in the case when at r = we have r,,,; = ^r°, that is the case when spatial 
positions are bounded in the same spatial scale in which their group velocities are bounded. 
In this case the evolution of the positions is described by the formula 



The rectilinear motion of positions of particle-like wavepackets is a direct consequence of the 
spatial homogeneity of the master system fll.ip . If the system were not spatially homogeneous, 
the motion of the positions of particle-like wavepackets would not be uniform, but we don't 
study that problem in this paper. In the rescaled coordinates y = gr the trajectory of every 
particle is a flxed, uniquely deflned straight line deflned uniquely if | — > as f), /5 — > 0. Notice 
that under above mentioned genericity condition the uniform and independent motion (11. 5p 
of the positions of all involved particle-like wavepackets {h k^,;, r,,;; r)} persists though 
they can collide in the space. In the latter case they simply pass through each other without 
signiflcant nonlinear interactions, and the nonlinear evolution with a high accuracy is reduced 
just to a nonlinear evolution of shapes of the particle-like wavepackets. In the case when the 
set of the principal wavevectors {k^;} satisfy certain resonance conditions some components 
of the original multi-particle wavepacket can evolve into a more complex structure which can 
be only partly localized in the space and, for instance, can be needle or pancake like. We do 
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not study in detail those more complex structures here. 

Now let us discuss in more detail the superposition principle introduced and studied for 
general multi-wavepackets in [9] in the particular case when initially all r*; = 0. Here we 
consider multi-particle wavepackets with arbitrary r^.^ and develop a new argument based on 
the analysis of an averaged wavepacket interaction system introduced in [8] . Assume that the 
initial data h for the evolution equation fll.ip to be the sum of a finite number of wavepackets 
(particle-like wavepackets) hi, I = 1, . . . , N, i.e. 

h = hi + ... + h^ (1.6) 

where the monochromaticity of every wavepacket is characterized by another small pa- 
rameter (3. The well known superposition principle is a fundamental property of every linear 
evolutionary system, stating that the solution U corresponding to the initial data h as in 
(11. 6p equals 

U = Ui + . . . + Ujv, for h = hi + . . . + hjv, (1.7) 

where U/ is the solution to the same linear problem with the initial data h/. 

Evidently the standard superposition principle can not hold exactly as a general principle 
for a nonlinear system, and, at the first glance, there is no expectation for it to hold even 
approximately. We show though that, in fact, the superposition principle does hold with a 
high accuracy for general dispersive nonlinear wave systems such as (II. ip provided that the 
initial data are a sum of generic particle-like wavepackets, and this constitutes one of the 
subjects of this paper. Namely, the superposition principle for nonlinear wave systems states 
that the solution U corresponding to the multi-particle wavepacket initial data h as in (II. 6p 
satisfies 

U = Ui + . . . + Uat + D, for h = hi + . . . + h^, where D is small. 

More detailed statement of the superposition principle for nonlinear evolution of wavepackets 
is as follows. We study the nonlinear evolution equation (II. ip on a finite time interval 

< r < r^, where > is a fixed number (1-8) 

which may depend on the L°° norm of the initial data h but, importantly, does not depend 
on g. We consider classes of initial data such that wave evolution governed by U.l\) is 
significantly nonlinear on time interval [0, r^,] and the effect of the nonlinearity F (U) does 
not vanish as g ^ 0. We assume that (3, g satisfy 

0< /3<1, 0<^<1, — <Ci with some Ci > 0. (1.9) 

The above condition of boundedness on the dispersion parameter ^ ensures that the dis- 
persive effects are not dominant and they do not suppress nonhnear effects, see [8], [9] for a 
discussion. 
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Let us introduce the solution operator S (h) (r) : h — U (r) relating the initial data h of 
the nonlinear evolution equation (II. ip to its solution U (t). Suppose that the initial state is 

a system of particle- like wavepackets or multi-particle wavepacket, namely h = ^^h^, with 
hi, I = 1, . . . , N being "generic" wavepackets. Then for all times < t < the following 
superposition principle holds 

^ (EL (^) = ^ (r) + D (r) , (1.10) 

II D (r)ll^ = sup II D (r)||^^ < Cs-f-r for any small 5 > 0. (1.11) 

0<r<r. P 

Obviously, the right-hand side of ( 11.111) may be small only if g < Cij3. There are ex- 
amples (see ^) in which D (r) is not small for q = Ci(3. In what follows we refer to a 
linear combination of particle-like wavepackets as a multi-particle wavepacket, and to sin- 
gle particle-like wavepackets which constitutes the multi-particle wavepacket as component 
particle wavepackets. 

Very often in theoretical studies of equations of the form (II. ip or ones reducible to it a 
functional dependence between q and j3 is imposed, resulting in a single small parameter. 
The most common scaling is ^ = The nonlinear evolution of wavepackets for a variety 
of equations which can be reduced to the form (II. ip was studied in numerous physical and 
mathematical papers, mostly by asymptotic expansions of solutions with respect to a single 
small parameter similar to /3, see |I3], [16], [20], [22], [27], [33], [31], [H], [15], [S], [H] 
and references therein. Often the asymptotic expansions are based on a specific ansatz 
prescribing a certain form to the solution. In our studies here we do not use asymptotic 
expansions with respect to a small parameter and do not prescribe a specific form to the 
solution, but we impose conditions on the initial data requiring it to be a wavepacket or a 
linear combination of wavepackets. Since we want to establish a general property of a wide 
class of systems, we apply a general enough dynamical approach. There is a number of general 
approaches developed for the studies of high- dimensional and infinite-dimensional nonlinear 
evolutionary systems of hyperbolic type, [I2], [15], [21], [2l], [3l], [37], [12], [17], [19], [5l], 
|53j ) and references therein. The approach we develop here is based on the introduction of 
a wavepacket interaction system. We show in [9] and here that solutions to this system are 
in a close relation to solutions of the original system. 

The superposition principle implies, in particular, that in the process of nonlinear evolu- 
tion every single wavepacket propagates almost independently of other wavepackets (even 
though they may "collide" in physical space for a certain period of time) and the exact 
solution equals the sum of particular single wavepacket solutions with a high precision. In 
particular, the dynamics of a solution with multi-wavepacket initial data is reduced to dy- 
namics of separate solutions with single wavepacket data. Note that the nonlinear evolution 
of a single wavepacket solution for many problems is studied in detail, namely it is well ap- 
proximated by its own nonlinear Schrodinger equation (NLS), see [20], [21], [31], [35], [17], 
|48j . |49j . [8] and references therein. 

Let us give now an elementary physical argument justifying the superposition principle 
which goes as follows. If there would be no nonlinearity, the system would be linear and. 
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consequently, the superposition principle would hold exactly. Hence, any deviation from it is 
due to the nonlinear interactions between wavepackets, and one has to estimate their impact. 
Suppose that initially at time r = the spatial extension s of every involved wavepacket 
is characterized by the parameter /3~^ as in (11.31) . Assume also (and it is quite an assump- 
tion) that the involved wavepackets evolving nonlinearly maintain somehow their wavepacket 
identities, including the group velocities and the spatial extensions. Then, consequently, the 
spatial extension of every involved wavepacket is propositional to and its group veloc- 
ity vi is proportional to The difference Av between any two different group velocities 
is also proportional to g^^. Then the time when two different wavepackets overlap in the 
space is proportional to s/ \ Av\ and, hence, to g/f3. Since the nonlinear term is of order one, 
the magnitude of the impact of the nonlinearity during this time interval should be roughly 
proportional to g/f3, which results in the same order of the magnitude of D in (I1.10l) - (ll.lll) . 
Observe, that this estimate is in agreement with our rigorous estimate of the magnitude of 
D in ffLTTD if we set there 5 = 0. 

The rigorous proof of the superposition principle presented here is not directly based 
on the above argument since it already implicitly relies on the principle. Though some 
components of the physical argument can be found in our rigorous proof. For example, 
we prove that the involved wavepackets maintain under the nonlinear evolution constant 
values of their wavevectors with well defined group velocities (the wavepacket preservation). 
The Theorem 16.121 allows to estimate spatial extensions of particle-like wavepackets under 
the nonlinear evolution. The proof of the superposition principle for general wavepackets 
provided in p] is based on general algebraic-functional considerations and on the theory 
of analytic operator expansions in Banach spaces. Here we develop an alternative approach 
with a proof based on properties of the wavepacket interaction systems introduced in [Sj . 

To provide a flexibility in formulating more specific statements related to the spatial 
localization of wavepackets we introduce a few types of wavepackets: 

• a single particle-like wavepacket w which is characterized by the following properties: (a) 
its modal decomposition involves only wavevectors from /3- vicinity of a single wavevector 
k^,, where /? > is a small parameter; (b) it is spatially localized in all directions and 
can be assigned its position r*; 

• a multi-particle wavepacket which is a system {wi} of particle-like wavepackets with 
the corresponding sets of wavevectors {k*^} and positions {r,,,^}; 

• a spatially localized multi-wavepacket which is a system {wi} with wi being either a 
particle- like wavepacket or a general wavepacket. 

We would like to note that a more detailed analysis, which is left for another paper, indi- 
cates that under certain resonance conditions nonlinear interactions of particle-like wavepack- 
ets may produce a spatially localized wavepacket w characterized by the following properties: 
(i) its modal decomposition involves only wavevectors from /3-vicinity of a single wavevector 
k^, where /3 > is a small parameter; (ii) it is only partly spatially localized in some, not 
necessarily all directions, and, for instance, it can be needle-like or pancake-like. 
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We also would like to point out that the particular form f 1 1.1 1) of dependence on the 
small parameter q is chosen so that appreciable nonlinear effects occur at times of order one. 
In fact, many important classes of problems involving small or parameters can be readily 
reduced to the framework of (11. II) by a simple rescaling. It can be seen from the following 
examples. The first example is a system with a small nonlinearity 

dtv = -iLv + af (v) , vj^^^ = h, < a < 1, (1.12) 

where the initial data is bounded uniformly in a. Such problems are reduced to (II. ip by 
the time rescaling r = ta. Note that here g = a and the finite time interval < r < 
corresponds to the long time interval <t < T^/a. 

The second example is a system with small initial data considered on long time intervals. 
The system itself has no small parameters but the initial data are small, namely 

dtv = —iLv + fo (v) , v|^^Q = aoh, < ao 1, where (1-13) 
fo (v) = (v) + f^-^^^ (v) + . . . , 

where ao is a small parameter and f-*"^ (v) is a homogeneous polynomial of degree m > 2. 
After rescaling v = aoV we obtain the following equation with a small nonlinearity 



dtV = -iLV + ar^i 



•^0 



f^) (V) + aof°(™+'^ (V) + . . . , V|,^o = h, (1.14) 



which is of the form of (11.121) with a = a™ ^. Introducing the slow time variable r = ta^ ^ 
we get from the above an equation of the form (II. ip . namely 

drV = -^LV + [f (V) + aof^™+'^ (V) + . . .] , = h, (1.15) 

where the nonlinearity does not vanish as — > 0. In this case g = and the finite time 

interval < r < corresponds to the long time interval < t < with small -C 1- 

The third example is related to a high-frequency carrier wave in the initial data. To be 
concrete, we consider the Nonlinear Schrodinger equation 

d^U = -idlU + ia \U\^ U, Ul^Q = hi {Mpx) e'^^^'^'' + /12 {Mpx) e^^^^*^^ + c.c, (1.16) 

where c.c. stands for complex conjugate of the prior term, and M ^ 1 is a large parameter. 
The equation (I1.16P can be readily recast into the form (ll.ip by change of variables y = Mr 
yielding 



drU = -i-dlU + la |[/|^ [/, = hi (pr) e'^*'' + /i2 (/9r) e''^*^" + c.c, (1.17) 



where g=^<^l. 



Summarizing the above analysis we list below important ingredients of our approach. 
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• The wave nonlinear evolution is analyzed based on the modal decomposition with re- 
spect to the linear part of the system. The significance of the modal decomposition 
to the nonlinear analysis is based on the following properties: (i) the wave modal 
amplitudes do not evolve under the linear evolution; (ii) the same amplitudes evolve 
slowly under the nonlinear evolution; (iii) modal decomposition is instrumental to the 
wavepacket definition including its spatial extension and the group velocity. 

• Components of multi-particle wavepacket are characterized by their wavevectors k^;, 
band numbers rii and spatial positions r,,;. The nonlinear evolution preserves k^,; and 
rii whereas the spatial positions evolve uniformly with the velocities ^Vu;„, (k*;). 

• The problem involves two small parameters (3 and g respectively in the initial data and 
coefficients of the master equation (11. ip . These parameters scale respectively (i) the 
range of wavevectors involved in its modal composition, with /3~^ scaling its spatial 
extension, and (ii) g scaling the ratio of the slow and the fast time scales. We make 
no assumption on the functional dependence between j3 and g, which are essentially 
independent and are subject only to inequalities. 

• The nonlinear evolution is studied for a finite time which may depend on, say, 
the amplitude of the initial excitation, and, importantly, r* is long enough to observe 
appreciable nonlinear phenomena which are not vanishingly small. The superposition 
principle can be extended to longer time intervals up to blow-up time or even infinity if 
relevant uniform in (3 and g estimates of solutions in appropriate norms are available. 

• In the chosen slow time scale there are two fast wave processes with typical time scale 
of order g which can be attributed to the linear operator L: (i) fast time oscillations 
resulting in time averaging and consequent suppression of many nonlinear interactions; 
(ii) fast wavepacket propagation with large group velocities resulting in effective weak- 
ening of nonlinear interactions which are not time-averaged because of resonances. It 
is these two processes provide mechanisms leading to the superposition principle. 

The rest of the paper is organized as follows. In the following Subsection 2.1 we introduce 
definitions of wavepackets, multiwavepackets and particle wavepackets. In subsection 2.1 we 
also formulate and briefiy discuss some important results of [8j which are used in this paper, 
and in Subsection 2.2 we formulate new results. In Section 3 we formulate conditions imposed 
on the linear and the nonlinear parts of the evolution equation (11.11) . and also introduce 
relevant concepts describing resonance interactions inside of wavepackets. In Section 4 we 
introduce an integral form of the basic evolution equation and study basic properties of 
involved operators. In Section 5 we introduce wavepacket interaction system describing 
the dynamics of wavepackets. In Section 6 we, first, define averaged wavepacket interaction 
system which plays a fundamental role in the analysis of the dynamics of multiwavepackets 
and then prove that solutions to this system approximate solutions to the original equation 
with high accuracy. We also discuss there properties of averaged nonlinearities, in particular, 
for universally and conditionally universal invariant wavepackets, and prove the fundamental 
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theorems on preservation of multi-particle wavepackets, namely Theorems 16. 131 and Theorem 
I2.1UI In Section 7 we prove the superposition principle using an approximate decoupling of 
the averaged wavepacket interaction system. In the last subsection of this section we prove 
some generalizations to the cases involving non-generic resonance interactions such as the 
second-harmonic and the third-harmonic generations. 

2 Statement of results 

This section consists of two subsections. In the first one we introduce basic concepts and 
terminology and formulate relevant results from [8] which are used latter on, and in the 
second one we formulate new results of this paper. 

2.1 Wavepackets and their basic properties 

Since the both linear operator L (— iV) and the nonlinearity F (U) are translation invariant, 
it is natural and convenient to recast the evolution equation f 1 1.1 1) by applying to it the Fourier 
transform with respect to the space variables r, namely 



and F is the Fourier form of the nonlinear operator F (U) involving convolutions, see (13. 9p 
for details. The equation (12.11) is written in terms of Fourier modes, and we call it the modal 
form of the original equation (11.11) . The most of our studies are conducted first for the modal 
form (12. ip of the evolution equation and carried over then to the original equation (II. ID . 

The nonlinear evolution equations (II. ip . (12. ip are commonly interpreted as describing 
wave propagation in a nonlinear medium. We assume that the linear part L (k) is a 2 J x 2 J 
Hermitian matrix with eigenvalues uJn,( (k) and eigenvectors g„^^ (k) satisfying 



L (k) g„,^ (k) = Un,c (k) g„,c (k) , C = ±, ^n,+ (k) > 0, Un,^ (k) < 0, n = 1, . . . , J, (2.3) 



where (k) are real-valued, continuous for all non-singular k functions, and vectors 
gn,c (k) G C^'^ have unit length in the standard Euclidean norm. The functions ct;„ ,^(k), 
n = 1, . . . , J, are called dispersion relations between the frequency u and the wavevector k 
with n being the band number. We assume that the eigenvalues are naturally ordered by 




(2.1) 



where U (k) is the Fourier transform of U (r), i.e. 




where r, k G M'^, 



(2.2) 



(k) > . . . > (k) > > (k) > . . . > cuj,_ (k) , 



(2.4) 
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and for almost every k (with respect to the standard Lebesgue measure) the eigenvalues are 
distinct and, consequently, the above inequalities become strict. Importantly, we also assume 
the following diagonal symmetry condition 

uJn-c (-k) = -ujn,(: (k) , C = ±, n = l,...,J, (2.5) 

which is naturally present in many physical problems (see also Remark 13.31 below), and is a 
fundamental condition imposed on the matrix L (k). Very often we use the abbreviation 

(k) = (k) . (2.6) 

In particular we obtain from 02.51) 

Un,^ (k) = -Un (-k) , UJn,i: (k) = (uJn (C^) , C = ±- (2-7) 

In addition to that in many examples we also have 

g„_^ (k) = _^ (^k) , where z* is complex conjugate to z. (2.8) 

We also use rather often the orthogonal projection n„^^ (k) in C^"' onto the complex line 
defined by the eigenvector g„^^ (k), namely 

n„,^ (k) u (k) = (k) gn,c (k) = un,c (k) , n = 1, . . . , J, C = ±. (2.9) 

As it is indicated by the title of this paper we study the nonlinear problem (II. ip for initial 
data h in the form of a properly defined particle-like wavepackets or, more generally, a sum 
of such wavepackets to which we refer as multi-particle wavepacket. The simplest example of 
a wavepacket w is provided by the following formula 

w (/3; r) = $+ il3 (r - r,)) e** )g„,+ (k,) , r G (2.10) 

where k* G M'^ is wavepacket principal wavevector, n is hand number, and /3 > is a small 
parameter. We refer to the pair (n, k,,,) in fl2.10p as wavepacket nk-pair and r* as wavepacket 
position. Observe that the space extension of the wavepacket w (/3; r) is proportional to (3~^ 
and it is large for small (3. Notice also that if /5 — > the wavepacket w (/3; r) as in fl2.10p 
tends, up to a constant factor, to the elementary eigenmode e'''* '"g„,^ (k*) of the operator 
L (— iV) with the corresponding eigenvalue Un,c_ (k*). We refer to wavepackets of the simple 
form fl2.10p as simple wavepackets to underline the very special way the parameter (3 enters its 
representation. The function $^ (r), which we call wavepacket envelope, describes its shape 
and it can be any scalar complex-valued regular enough function, for example a function 
from Schwartz space. Importantly, as /? — ^ the L°° norm of a wavepacket ^KJW remains 
constant, hence nonlinear effects in U.l\) remain strong. 

Evolution of wavepackets in problems which can be reduced to the form (11. ip were studied 
for a variety of equations in numerous physical and mathematical papers, mostly by asymp- 
totic expansions with respect to a single small parameter similar to /3, see [13], p^, [20j . 
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[22], [2Z1, [23], [M], [S], [15], [IZ], [IB] and references therein. We are interested in general 
properties of evolutionary systems of the form f 1 1.11) with wavepacket initial data which hold 
for a wide class of nonlinearities and all values of the space dimensions d and the number 2 J 
of the system components. Our approach is not based on asymptotic expansions but involves 
the two small parameters (3 and q with mild constraints (11.91) on their relative smallness. 
The constraints can be expressed either in the form of certain inequalities or equalities, and 
a possible simple form of such a constraint can be a power law 

j3 = Cg^ where C > and x > are arbitrary constants. (2-11) 

Of course, general features of wavepacket evolution are independent of particular values of the 
constant C. In addition to that, some fundamental properties such as wavepacket preservation 
are also totally independent on particular choice of the values of x in (12. lip , whereas other 
properties are independent of x as it varies in certain intervals. For for instance, dispersion 
effects are dominant for x < 1/2, whereas the wavepacket superposition principle of [8] holds 
for X < 1. 

To eliminate unbounded (as ^ — > 0) linear term in (12. ip by replacing it with a highly 
oscillatory factor we introduce the slow variable u (k, r) by the formula 

U (k, r) = e-^^^'^^u (k, r) , (2.12) 
and get the following equation for u (k, r) 

drVi = ef^F (^e^'^uj , u|^^o = h, (2.13) 

which, in turn, can be transformed by time integration into the integral form 

u = J- (u) + h, J- (u) = ^ e^^F (^e^^u (r')) dr' (2.14) 

with explicitly defined nonlinear polynomial integral operator JF = JF(^). This operator is 
bounded uniformly with respect to g in the Banach space E = C {[0, r*] , L^) .This space has 
functions v (k, r), < r < as elements and has the norm 

||v(k,r)||^=||v(k,r)||^(p,.,],Li)= sup / |v(k,r)|dk, (2.15) 

0<T<T, JlRd 

where is the Lebesgue space of functions v (k) with the standard norm 

||v(.)ILi= / |v(k)|dk. (2.16) 

JR'' 

Sometimes we use more general weighted spaces L^'" with the norm 

= / (l + |k|)'^|v(k)| dk, a>0. (2.17) 
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The space C ([0, r*] , L^''') with the norm 

||v(k,r)||^^= sup / (l + |k|r|v(k,r)|dk. (2.18) 

0<r<r* J^d 

is denoted by Ea, and, obviously, Eq = E. 

A rather elementary existence and uniqueness theorem (Theorem 14.81) implies that if 
h G L^'" then for a small and, importantly, independent of g constant r* > this equation 
has a unique solution 

n{T)=g (g) , h) (r) , r G [0, r,] , u G C ([0, r,] , L'''^) , (2.19) 

where Q denotes the solution operator for the equation (12.141) . If u (k, r) is a solution to 
the equation (I2.14p we call the function U (r, r) defined by (I2.12p . (12. 2p an F-solution to 
the equation (II. ip . We denote by V" the space of functions V (r) such that their Fourier 

transform V (k) belongs to L^, and define ||V|||^i = V . Since 

||V||^^ < i^Tiy^ V and C L°°, (2.20) 



F-solutions to belong to (^[0, r,] , L^) C ([0, r,] , L°°). 

We would like to define wavepackets in a form which explicitly allows them to be real 
valued. This is accomplished based on the symmetry (12.50 of the dispersion relations, which 
allows to introduce a doublet wavepacket 

w (/5; r) = $+ iP (r - r,)) e^'^*-('-'*)g„,+ (k,) + $_ (/5 (r - r,)) e"'''* )g„,_ (-k,) . (2.21) 

Such a wavepacket is real if $_ (r), g„ _. (— k^,) are complex conjugate respectively to (r), 
g„,+ (k*), i.e. if 

$_ (r) = $; (r) , g„,+ (k,) = g„,_ (-k,)* . (2.22) 

Usually when considering wavepackets with n/c-pair {n, k^,) we mean doublet ones as in 
(I2.2ip . but sometimes we use the term wavepacket also for an elementary one as defined by 
(I2.10p . Note that the latter use is consistent with the former one since it is possible to take 
one of two terms in (I2.2ip to be zero. 

Below we give a precise definition of a wavepacket. To identify characteristic properties 
of a wavepacket suitable for our needs, let us look at the Fourier transform w (/5; k) of an 
elementary wavepacket w (/5; r) defined by (I2.10p . that is 

w (/3; k) = /3-'^e"* '-*<l {p-' (k - k,)) g„,c (k,) . (2.23) 

We call such w(/?;k) wavepacket too and it possesses the following properties: (i) its 
norm is bounded (in fact, constant), uniformly in /3 — >■ 0; (ii) for every e > the value 
w(/5;k) for every k outside a /?^~'^-neighborhood of k*, and the convergence is faster 
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than any power of /5 if $ is a Schwartz function. To exphcitly interpret the last property we 
introduce a cutoff function \E' {rj) which is infinitely smooth and such that 

^ (rj) > 0, ^ (t;) = 1 for \r)\ < 1/2, $ (rj) = for |t7| > 1, (2.24) 

and its shifted/rescaled modification 

^ k,; k) = ^ (k - k,)) . (2.25) 

If an elementary wavepacket w (/5; r) is defined by (I2.23P with $ (r) being a Schwartz function 
then 

II (1 - VI/ (/5I-, k,; ■)) w if3; Oil < < /? < 1, (2.26) 

and the inequality holds for arbitrarily small e > and arbitrarily large s > 0. Based on the 
above discussion we give the following definition of a wavepacket which is a minor variation 
of [HI Definiton 8]. 

Definition 2.1 (single-band wavepacket) Let e be a fixed number, < e < 1. For a 
given band number n G {1, . . . , J} and a principal wavevector k* e M*^ a function h (/5; k) 
is called a wavepacket with rafc-pair (n, k*) and the degree of regularity s > if for small 
(3 < with some (3^ > it satisfies the following conditions: (i) h(/?;k) is L^-bounded 
uniformly in [3, i.e. 

h (/?; ■) < C". < /3 < /3o for some C > 0; (2.27) 

(ii) h (/3; k) is composed essentially of two functions (/3; k), ^ = ±, which take values in 
n-th band eigenspace ofL (k) and are localized near ^k*, namely 

h {13- k) = h_ [13- k) + h+ {I3- k) + D,, < /5 < 13,, (2.28) 

where the components h± (/5; k) satisfy the condition 

\ (/5; k) = XI/ (/3I-72, Ck,; k) n„,c (k) (/3; k) , C = ±, (2.29) 

where ^ (■, Ck*, /?^~^) zs defined by h2.25\] and Dh is small, namely it satisfies the inequality 

WDhWii < C'I3\ 0<p <Pq, for some C > 0. (2.30) 

The inverse Fourier transform h (/5; r) of a wavepacket h (/9; k) is also called a wavepacket. 

Evidently, if a wavepacket has the degree of regularity s, it also has a smaller degree of 
regularity s' < s with the same e. Observe that the degree of regularity s is related to the 
smoothness of (r) as in (12.101) so that the higher is the smoothness the higher | can be 
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taken. Namely, if $^ G L}'"" then one can take in (I2.30p any s < ae according to the following 
inequahty: 



dr)<(3'' 



< C(3'. 



;i-vi>(/5^T7))$c ii) 

For example, if we define similarly to (12.291) and fl2.23p by the formula 
he k) = VP (k - k,)) (k - k,)) n„,e (k) g 



(2.3i; 



(2.32) 



where $e (k) is a scalar Schwartz function and g is a vector, then according to( l2.3ip estimate 
(I2.30p holds and {f3\ k) is a wavepacket with arbitrarily large degree of regularity s for 
any given e such that < e < 1. 

Now let us define a particle-like wavepacket following to the ideas indicated in the Intro- 
duction. 

Definition 2.2 ( single-band particle-like wavepacket) We call a function h(/3;k) = 
h (/5, r*; k), G M!^ , a particle-like wavepacket with the position r*, nk-pair (ra, k,,,) and the 
degree of regularity s > if (i) for every it is a wavepacket with the degree of regularity s 
in the sense of the above Definition \2. 1\ with constants C,C' independent ofr^ G M.'^; (ii) 
in 112. 2S\) satisfy the inequalities 



Vk(e''-*^haAr,;k)) dk < C,P 



(2.33) 



where Ci > is an independent of (3 and constant, e is the same as in Definition \2.1\ 
The inverse Fourier transform h (/?; r) of a wavepacket h (/5; k) is also called a particle-like 
wavepacket with the position r^. We also introduce quantity 



(2.34) 



which we refer to as the position detection function for the wavepacket h (/5, r*; k). 
Note that the left-hand side of (12.330 coincides with a ( r^,, (r^,) j . 



Remark 2.3 // h (/5; k) = h(/9, r^,;k) is a particle-like wavepacket with a position r* then, 
applying inverse Fourier transform to (/?, r*; k) and Vkh^ (/5, r^.; k) as in Ii2.2\) we obtain 
a function h (/5, r*; r) which satisfies 

(2.35) 



|r - r*| |h(; (/?, r*;r)| < (27r) a(^r*,h, 

implying that (/9;i")| < a ^r^,,hej |r — r^,!""*^. This inequality is useful for large |r — r*|, 
whereas for bounded |r — r*| l[2.21^ implies a simpler inequality 



W (/?,r.;r)| < (27r)- 
14 



< C. 



(2.36) 



Inequalities \2. 35\) and \2. 33\) suggest that the quantity a (^r,,,, (r^,) j can be interpreted as 
a size of the particle-like wavepacket r^.; k). 

Evidently a particle-like wavepacket is a wave and not a point. Hence the above definition 
of its position has a degree of uncertainty, allowing, for example, to replace r* by r,, + a with a 
fixed vector a (but not allowing unbounded values of a). The above definition of particle-like 
wavepacket position was crafted to meet the following requirements: (i) a system of particle- 
like wavepackets remains to be such a system under the nonlinear evolution; (ii) it is possible 
(in an appropriate scale) to describe the trajectories traced out by the positions of a system 
of particle-wavepackets 

Remark 2.4 Typical dependence of the inverse Fourier transform h (/3, r,,; r) of a wavepacket 
h (/3, r*; k) on r* is provided by spatial shifts by r^: as in h2.21\) . namely 

h(/3,r,;r) = $(/3 (r - r,)) e^'^-^^-^'-^g 

with a constant g. For such a function h and for any r'^ G M°' 



a (<,h(r. 



Vk(r"e''^"-'*h(Ar,;k) 
g 



d^ikr' — ikr. 



<l(k) 



|g| 



Hence, taking for simplicity 

1 



i(r:-r,)$(k') + -^Vk'$(k') 

= 1, we obtain 









1 






)>- 




$ 


Li ~ ^ 







(2.37) 



For small |r'^ — r*| ^ ^ we see that the position detection function a ^r'^,hj is of order 
O which is in the agreement with \2.3'J^) . For large |r'^ — r^,| ^ ^ the a ^r'^,hj is 

approximately proportional to |r'^ — r*|. Therefore if we know a ^r'^, h {r^)^ as a function of 

v'^ we can recover the value of r^, with the accuracy of order O (z?""^'*^) with arbitrary small 
e. Namely, let us take arbitrary small e > and some C > and consider the set 



: a\ r'^, h (r* 



< 



c 



(2.38) 



which should provide an approximate location ofr^. According to \2.31^ , lies in this set 
for small (3. If r'^ lies in this set then 



C/3- 



-l-e 



> a ( r' hfr* 









1 






)>- 










Li 
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and \r'^ — < Ci/3 ^ + C2P ■ Hence the diameter of the B {(3) is of order O {^(3 ^ . 
Observe, taking into account Remark \2.'J\ that the accuracy of the wavepacket location ob- 
viously cannot be better than its size a ^r*,hf (r*)j ~ P^'^ . The above analysis suggest that 

the function h r^,; r) can be viewed as pseudoshifts of the function h 0; r) by vectors 
r* G M'^ in the sense that the regular spatial shift by r* is combined with a variation of the 
shape o/h(/3,0;r) which is limited by the fundamental condition Ii2.33\) . In other words, 
according Definition \2.2\ as wavepacecket moves from to by a corresponding spatial shift 
it is allowed to change its shape subject to the fundamental condition Ii2.33\) . The later is 
instrumental for capturing nonlinear evolution of particle-like wavepackets governed by an 
equation of the form U.l\) . 

Remark 2.5 The set B {(3) defined by Ii2. 38\) gives an approximate location of the support of 
the function h {j3, r^,; k) not only in the special case considered in Remark \2.4\ but also when 
h r^,; r) is a general particle-like wavepacket. One can apply with obvious modifications the 
above argument for e^^^*h {/3, r*; k) in place of $ (k) using h2.33\) . Here we give an alternative 

argument based on l{2.35\) . Notice that condition a ^r^,o,h(r*)j < 0(3^^^^ obviously can be 

satisfied not only by r^,o = r^,. But one can show that the diameter of the set of such r^,o is 
estimated by O (^P^^^'^Y Indeed, assume that a given function h(/?,r) does not vanish at a 
given point tq, that is |h ro)| > Cq > for all [3 < Pq. The fulfillment of 112. 33\) for the 
function h {P, r) with two different values of r^, namely = r'^ and = r" implies that 

a (r:,h) < Cl/3-^-^ a (<, h) < C2p-'-\ 

and according to h2. 35^) for all r 

Ir-r'J |h(/3,r)| < {2'n)-'' C,p-^-\ |r - |h(Ar)| < {27,)-" C^P'^-^ 

Hence, 

ko - rj < , |ro - rj < , 



and 



Co Co 



V.-<\<c^P' 



Note that if we rescale variables r and r* as in Example 12. 131 namely ^?r = y and fpr* = y,,, 
with Q = P'^, the diameter of the set B (P) in y- coordinates is of order ^ 1, and, 
hence, this set gives a good approximation for the location of the particle-like wavepacket as 
/3 — s> 0. It is important to notice, that our method to locate the support of wavepackets is 
applicable to very general wavepackets and does not use their specific form. This flexibility 
allows us to prove that particle-like wavepackets and their positions are well defined during 
nonlinear dynamics of generic equations with rather general initial data which form infinite- 
dimensional function spaces. Another approaches to describe dynamics of waves are applied 
to situations where solutions under considerations can be parametrized by a finite number of 
parameters and the dynamics of the parameters describes dynamics of the solutions. See for 
example [25], [26] where dynamics of centers of solitions is described. 
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Remark 2.6 Note that for a single wavepacket initial data h r — r'^) one can make a 
change of variables to a moving frame (x, r), namely (r, r) = (x + vr, r), where v = 
^Vcu (k,,) is the group velocity; this change of variables makes the group velocity zero. Often 
it is possible to prove that dynamics preserves functions which decay at infinity , namely if 
the initial data h {(3, x) decays at the spatial infinity then the solution U (/?, x, r) also decays 
at infinity (though corresponding proofs can be rather technical). This property can be re- 
formulated in rescaled y variables as follows: if initial data are localized about zero, then the 
solution is localized about zero as well. Then, using the fact that the equation has constant 
coefficients, we observe that the solution U y — y'^, r), corresponding to h(/5, y — y^), is 
localized about y'^ provided that h y) was localized about the origin. Note that in this 
paper we consider much more complicated case of multiple wavepackets. Even in the sim- 
plest case of the initial multiwavepacket which involves only two components, namely the 
wavepacket h r) = hi r — r'J + h2 r — r'^') with two principal wave vectors ki* ^ k2* 
one evidently cannot use the above considerations based on the change of variables and the 
translational invarance. Using other arguments devoloped in this paper we prove that systems 
of particle-like wavepackets remain localized in the process of the nonlinear evolution. 

Note that similarly to fll.2l) and fll.4p a function of the form 

k-k 



h 



defined for any pair of r*i and r*2 where h is a Schwartz function and all constants in 
Definition 12.11 are independent of r,.i,r^,2 G M.^, is not a single particle-like wavepacket since 
it does not have a single wavepacket position r,., but rather it is a sum of two particle-like 
wavepackets with two positions r,,i and r,,2. 

We want to emphasize once more that a particle-like wavepacket is defined as the family 
h r^; k) with being an independent variable running the entire space M'^, see, for ex- 
ample, (11.21) . (11. 3p and (I2.2ip . In particular, we can choose a dependence of r,, on (3 and g. 
An interesting type of such a dependence is r* = g where g satisfies (12.111) as we discuss 
below in the Example 12.131 

Our special interest is in the waves that are finite sums of wavepackets which we refer to 
as multi-wavepackets. 

Definition 2.7 (multi-wavepacket) Let S be a set of nk-pairs: 

S = {{ni, Ki) , / = 1, . . . , AT} c S = {1, . . . , J} X {m, Ki) ^ {nu, Ku) for I I', (2.39) 

and N = \S\ be their number. Let Ks be a set consisting of all different wavevectors k*; 
involved in S with \Ks\ < N being the number of its elements. Ks is called wavepacket 
/c-spectrum and without loss of genericity we assume the indexing of elements {ni, k^,;) in S 
to be such that 

Ks = {K^,^ = l,...AKs\}, ^.e. I = t for 1 < t < \Ks\ . (2.40) 
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A function h {(3) = h k) is called a multi- wavepacket with nfc-spectrum S if it is a finite 
sum of wavepackets, namely 



N 



h {(3; k) = ^ hi k) , < (3 < (3^ for some (3^ > 0, (2.41) 



1=1 



where hi, I = 1, . . . , N , is a wavepacket with nk-pair {k^,i,ni) E S as in Definition \2.1\ If 
all the wavepackets h; k) = h/ v^.i; k) are particle-like ones with respective positions r*; 
then the multi-wavepacket is called muhi-particle wavepacket and we refer to (r*i, . . . , r^,7v) 
as its position vector. 

Note that if h(/9;k) is a wavepacket then h(/?;k) + O {{3^) is also a wavepacket with 
the same n/c-spectrum, and the same is true for muhi-wavepackets. Hence, we can introduce 
multi- wavepackets equivalence relation "~" of the degree s by 

hi k) ~ ha k) if hi k) - ha k) < Cl3' for some constant C > 0. (2.42) 

Note that condition (12.331) does not impose restrictions on the term Dh in (I2.28p . therefore 
this equivalence can be applied to particle wavepackets. 

Let us turn now to the abstract nonlinear problem (12.141) where (i) JF = JF [q) depends on 
Q and (ii) the initial data h = h {(3) is a multi-wavepacket depending on (3. We would like to 
state our first theorem on multi-wavepacket preservation under the evolution (12.141) for (3, g —>■ 
0, which holds provided its n/c-spectrum S satisfies a natural condition called resonance 
invariance. This condition is intimately related to the so-called phase and frequency matching 
conditions for stronger nonlinear interactions, and its concise formulation is as follows. We 
define for given dispersion relations {un (k)} and any finite set S G {1, . . . , J} x M*^ another 
finite set n{S) C {1, . . . , J} x R'^ where 7^ is a certain algebraic operation described in 
Definition 13.81 below. It turns out that for any S always S* C 7^ (5*), but if IZ {S) = S we call 
S resonance invariant. The condition of resonance invariance is instrumental for the multi- 
wavepacket preservation, and there are examples showing that if it fails, i.e. TZ {S) ^ S, 
the wavepacket preservation does not hold. Importantly, the resonance invariance R (S) = 
S allows resonances inside the multi-wavepacket, that includes, in particular, resonances 
associated with the second and the third harmonic generations, resonant four- wave interaction 
etc. In this paper we will use basic results on wavepacket preservation obtained in [8] , 
and we formulate theorems from [8] we need here. Since we use constructions from [8], for 
completeness we provide also their proofs in the following subsections. The following two 
theorems are proved in [8j. 

Theorem 2.8 (multi-wavepacket preservation) Suppose that the nonlinear evolution is 
governed by ^2.14^ and the initial data h = h (/5; k) is a multi-wavepacket with nk-spectrum 
S and the regularity degree s, and assume S to he resonance invariant (see Definition \ 3.^ 
below). Let p{P) be any function satisfying 

0< p{l3) < Cl3', for some constant C > 0, (2.43) 
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and let us set q = p{(3). Then the solution u(r,/?) = Q |^JF (p (/?)), h (r) to ^2.14\ ) for 
any t G [0, r*] is a multi-wavepacket with nk-spectrum S and the regularity degree s, i.e. 

u (r, /3; k) = ^ (r, /5; k) , where u; is wavepacket with nk-pair {ni, k*;) G 5. (2.44) 

r/ie time interval length r^, > depends only on V- -norms of hi (/3; k) and N. The presen- 
tation ^2.44^ is unique up to the equivalence (2^4^ of degree s. 

The above statement can be interpreted as follows. Modes in n/c-spectrum S are always 
resonance coupled with modes in TZ {S) through the nonlinear interactions, but if TZ {S) = S 
then (i) all resonance interactions occur inside S and (ii) only small vicinity of 5* is involved 
in nonlinear interactions leading to the multi-wavepacket preservation. 

The statement of Theorems 12.81 directly follows from the following general theorem proved 
in [H]. 

Theorem 2.9 (multi-wavepacket approximation) Let the initial data h in the integral 
equation ^2. 14\ ) be a multi-wavepacket h (/5; k) with nk-spectrum S as in Ii2. 39\) . the regularity 
degree s and with the parameter e > as in Definition \2.1\ Assume that S is resonance 
invariant in the sense of Definition \3.8\ below. Let the cutoff function ^ (/5 k^,; k) and the 
eigenvector projectors Iln,± (k) be defined by l\2.25\) and l\2. 9\) respectively. For a solution u 
of i2.14\j we set 

{(3- r, k) = [5^^^^ v]/ Ck,,; k) n„,,^ (k)] u (/3; r, k) , Z = 1, . . . , iV. (2.45) 

Then every such (/5; r, k) is a wavepacket and 



sup 

0<T<T, 



u(/3;r,k)->;f>K/5;^,k) ^^<C^Q + C2(3' (2.46) 



where the constants C, Ci do not depend on e, s and (3, and the constant C2 does not depend 
on (3 and g. 

We would like to point out also that Theorem 12.81 allows to take values u (r*) as new 
wavepacket initial data for fll.ip and extend the wavepacket invariance of a solution to the 
next time interval < r < r^i. This observation allows to extend the wavepacket invariance 
to larger values of r (up to blow-up time or infinity) if some additional information about 
solutions with wavepacket initial data is available, see [8]. 

Note that the wavepacket form of solutions can be used to obtain long-time estimates 
of solutions. Namely, very often behavior of every single wavepacket is well approximated 
by its own nonlinear Schrodinger equation (NLS), see [12], [3H], ED], |2Z], [M], [SS], [IS], 
[17], [IH], [IS] and references therein, see also Section 6. Many features of the dynamics 
governed by NLS-type equations are well-understood, see [17], [18], [39], [16], [50], [51] and 
references therein. These results can be used to obtain long-time estimates for every single 
wavepacket (as, for example, in [35]) and, with the help of the superposition principle, for 
the multiwavepacket solution. 
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2.2 Formulation of new results on particle wavepackets 

In this paper we prove the following refinement of Theorem 12.81 for the case of multi-particle 
wavepackets. 

Theorem 2.10 (multi-particle wavepacket preservation) Assume that conditions of The- 
orem \2.9[ hold and, in addition to that, the initial data h = h (/5; k) is a multi-particle 
wavepacket of degree s with positions r^,i, . . . , r^^r and the multi-particle wavepacket is uni- 
versally resonance invariant in the sense of Definition \3. Assume also that 

p{f3)<C(3'\s,>{). (2.47) 



Then the solution u(/5;r) = Q ^JF (p (/5)) , h (/5) j (r) to \2.14\) for any r G [0, r 



IS a 



multi-particle wavepacket with the same nk-spectrum S and the same positions r*i, . . . , r^,7v- 
Namely, ^2.46\ ) holds where ui is wavepacket with nk-pair {ni,'k^i) G 5* defined by 11^2.45^ , 
the constants C,Ci,C2 do not depend on r^i. and every u; is equivalent in the sense of the 
equivalence \2.4^ of degree Si = min (s, Sq) to a particle wavepacket with the position r^[. 

Remark 2.11 Note that in the statement of the above theorem the positions r^,i, . . . , r^^r of 
wavepackets which compose the solution u (/5; r, k) of l{2.13\) and \2.14^ do not depend on r 
and, hence, do not move. Note also that the solution U (/?; r, k) of the original equation \2. 
related to u (/5; r, k) by the change of variables h2.12\} . is composed of wavepackets (/?; r, r), 
corresponding to Ui{f3;T,r), have their positions moving with respective constant velocities 
VfcCiJ (k^i) (see for details Remark 4-1, see also the following corollary). 



Using Proposition 14.21 we obtain from Theorem 12.101 the following corollary. 

Corollary 2.12 Let conditions of Theorem \2.1U\ hold and U (/?; r, k) be defined by l[2.12\) in 
terms of u{f3]T,k). Let 

11 

— < C, with some C, < /3 < 0< g< -. (2.48) 

Then U (/?; r, k) is for every r G [0,t^:] a particle multi-wavepacket in the sense of Defini- 
tion \2.2\ with the same nk-spectrum S , regularity si and with r-dependent positions r*/ + 

In the following example we consider the case where spatial positions of wavepackets have 
a specific dependence on parameter g, namely r^ = r^/ g. 

Example 2.13 (Wavepacket trajectories and collisions) Let us rescale the coordinates 
in the physical space as follows 

gr = y (2.49) 
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with the consequent resettling of the wttvevector vttrittble (duttl with respect to Fourier trttns- 
form) k = grj. It follows then thttt under the evolution U.l\) the group velocity of tt Wttvepacket 
with tt Wttvevector k* in the new coordintttes y becomes Vfco; (k^,) ttnd evidently is of order 
one. If we set the positions r^i = g with fixed r^i, then ttccording to h2. 5*5]) wttvepttckets 

|h(/9;r)| in y-vttrittbles httve chttrttcteristic sptttittl settle y — r*; ~ go. ^i"*«,hj ~ gP^^ which 

is smttll if g/P is smttll. The positions of pttrticle-like wttvepttckets (quttsipttrticles) U(y/f?, t) 
are initittlly locttted ttt yi = r°; ttnd propttgttte with the group velocities VA:Ci;(k*;). Their 
trajectories are strttight lines in the spttce M*^ described by 

y = rVkUJ (Ki) + rli, < r < 

(compttre with M.5\) ). The trttjectories mtty intersect, indicttting "collisions" of quttsipttrticles. 
Our results (Theorem \2.10\) show thttt if tt multi-pttrticle wttvepttcket initittlly wtts universttlly 
resonttnce invttrittnt, then the involved pttrticle-like wttvepttckets preserve their identity in spite 
of collisions ttnd the fttct thttt the nonlinettr interttctions with other wttvepttckets (quttsipttr- 
ticles) ttre not smttll, if fttct, they ttre of order one. Note thttt can be chosen ttrbitrttrily 
implying thttt up to N {N — 1) collisions cttn occur on the time intervttl [0, r*] on which we 
study the system evolution. 

To formulate the approximate superposition principle for multi-particle wavepackets, we 
introduce now the solution opertttor Q mapping the initial data h into the solution U = ^ 

of the modal evolution equation fl2.14p . This operator is defined for h < -R according to 

the existence and uniqueness Theorem 14. 7[ The main result of this paper is the following 
statement. 



Theorem 2.14 (superposition principle ) Suppose thttt the initittl dttttt h of \2.14^ 
multi-pttrticle wttvepttcket of the form 



IS tt 



N 

h = ^h;, iVmax h; < ^, (2.50) 



1=1 



stttisfying Definition \2. 7| ttnd its nk-spectrum is universttlly resonttnce invttrittnt in the sense 
of Definition \3.S\ Suppose ttlso thttt thttt the group velocities of wttvepttckets ttre different, 
nttmely 

WkUJni^ (Kh) 7^ Vk^^n,^ (Kh) ifh 7^ h (2.51) 



ttnd thttt I!i2.4^ holds. Then the solution u = Q {hj to the evolution equtttion I!i2.14\ ) stttisfies 
the following ttpproximttte superposition principle 

^ E^O W+D, (2.52) 

\i=i J 1=1 
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with a small remainder D (t) such that 



sup 

0<r<r, 



DH 



\\np\ 



(2.53) 



where (i) e is the same as in Definition \2. 1\ and can be arbitrary small; (ii) does not depend 
on P, Q, v^i and e; (Hi) does not depend on (3, g and positions r*; . 

A particular case of the above Theorem in which there was no dependence on r*; was 
proved in [9] by a different method based on the theory of analytic operators in Banach 
spaces. The condition (12.511) can be relaxed if the initial positions of involved particle-like 
wavepackets are far apart, and the the corresponding results are formulated in the theorem 
below and in Example 12. 131 



Theorem 2.15 (superposition principle ) Suppose that the initial data h of ^2.14\ ) 
a multi-particle wavepacket of the form ^2.50) with a universally resonance invariant nk- 
spectrum in the sense of Definition \3.8\ and ^2.48^ holds . Suppose also that either the group 
velocities of wavepackets are different, namely ^2.51\) holds. or the positions r^i satisfy the 
inequality 



I ^*/2 I 



< 



(2.54) 



where the constant C^^2 is the same as in ^3.2\) . Then the solution u = Q |^hj to the evolution 
equation {2.14) satisfies the approximate superposition principle ^2. 52\} . ^2.53^ . 



We prove in this paper further generalizations of the particle-like wavepacket preservation 
and the superposition principle to the cases where the nfc-spectrum of a multi-wavepacket 
is not universal resonance invariant such as the cases of multi-wavepackets involving the sec- 
ond and the third harmonic generation. In particular, we prove Theorem 17.51 showing that 
that many (but, may be, not all) components of involved wavepackets remain spatially local- 
ized. Another Theorem 17.71 extends the superposition principle to the case when resonance 
interactions between components of a multi-wavepackets can occur. 



3 Conditions and definitions 

In this section we formulate and discuss all definitions and conditions under which we study 
the nonlinear evolutionary system (11.11) through its modal, Fourier form (12. ip . Most of the 
conditions and definitions are naturally formulated for the modal form (12.11) . and this is one 
of the reasons we use it as the basic one. 
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3.1 Linear part 



The basic properties of the hnear part L (k) of the system (12.11) . which is a 2 J x 2 J Hermitian 
matrix with eigenvalues u;„,c (k) , has been already discussed in the Introduction. To account 
for all needed properties of L (k) we define the singular set of points k. 

Definition 3.1 (band-crossing points) We callkQ a band-crossing point /or L (k) ifuJn+i,(^O^o) = 
^n,c (ko) for some n,( orL (k) is not continuous at kg or if Ui^± (ko) = 0. The set of such 
points is denoted by Ubc- 

In the next Condition we collect all constraints imposed on the linear operator L (k). 

Condition 3.2 (linear part) The linear part L (k) of the system 1(2. 1\) is a 2 J x 2J Her- 
mitian matrix with eigenvalues ujn,c (k) and corresponding eigenvectors g„_^ (k) satisfying for 
k ^ (Tbc the basic relations li2.3\) - l[27S\) . In addition to that we assume: 

(i) the set of band-crossing points Ubc is a closed, nowhere dense set in and has zero 
Lebesgue measure; 

(a) the entries of the Hermitian matrix L (k) are infinitely difjerentiable in k for all k ^ Ubc 
that readily implies via the spectral theory, f^ . infinite differentiability of all eigenval- 
ues ujn (k) in k for all k ^ cTbc/ 

(Hi) L (k) satisfies a polynomial bound 

||L(k)|| <C{1 + |kn , k G for some C > andp> 0. (3.1) 



Note that since ujn,(; (k) are smooth if k ^ cTbc the following relations hold: 

max |VkC^„,,^| < max |Vkt^„,,J < C^,2, (3.2) 

|k±k,;|<7ro, «=l,...,Ar, |k±k,,|<7ro, «=l,.-,Af, 



where C^^^i and C^_2 are positive constants and 



TTo = min min(dist{±k*/,crbc},l)- (3.3) 

2 1=1,. ..,N 



Remark 3.3 (dispersion relations symmetry) The symmetry condition Ii2.5\) on the dis- 
persion relations naturally arise in many physical problems, for example Maxwell equations 
in periodic media, see JEl-l^, JE/, or when L (k) originates from a Hamiltonian. We would 
like to stress that this symmetry conditions are not imposed to simplify studies but rather to 
take into account fundamental symmetries of physical media. The symmetry causes resonant 
nonlinear interactions, which create non-trivial effects. Interestingly, many problems without 
symmetries can be put into the framework with the symmetry by a certain extension, f^. 
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Remark 3.4 (band-crossing points) Band-crossing points are discussed in more details 
in Section 5.4], JM Sections 4-i, 4-^]- In particular, generically the setai,^ of band- crossing 
point is a manifold of the dimension d — 2. Notice also that there is an natural ambiguity in 
the definition of a normalized eigenvector gn.c (k) of L (k) which is defined up to a complex 
number ^ with \^\ = 1. This ambiguity may not allow an eigenvector g„.^ (k) which can be 
a locally smooth function in k to be a uniquely defined continuous function in k globally for 
all k ^ (Tbc because of a possibility of branching. But, importantly, the orthogonal projector 
Iln,(^ (k) on g„^^ (k) as defined by l{2.9\) is uniquely defined and, consequently, infinitely dif- 
ferentiable in k via the spectral theory, fSB], for all k ^ cibc- Since we consider U (k) as an 
element of the space and cxbc is of zero Lebesgue measure considering k ^ cxbc is sufficient 
for us. 

We introduce for vectors u G C^"' their expansion with respect to the orthonormal 
basis{g„,^ (k)}: 

J .7 

^(^) = jm g„,c i^) = J2J2 , <C (k) = n„,c (k) u (k) (3.4) 

n=l C=± n=l C=± 

and we refer to it as the modal decomposition of u (k) and to ,^ (k) as the modal coefficients 
of u (k). Evidently 

y~^] ' ^ _|_ n„_^ (k) = I2J, wliere I2J is the 2 J x 2 J identity matrix. (3.5) 

Notice that we can define the action of the operator L (— iVr) on any Schwartz function Y (r) 
by the formula 

L(^^)Y(k) = L(k)Y(k), (3.6) 

where in view of the polynomial bound (13.11) the order of L does not exceed p. In a special 
case when all the entries of L (k) are polynomials (13.61) turns into the action of the differential 
operator with constant coefficients. 

3.2 Nonlinear part 

The nonlinear term F in (12. ip is assumed to be a general functional polynomial of the form 

F (tj) = ("U"^) , where F*-™"-* is m-homogeneous polylinear operator, (3.7) 

Tip = {mi, . . . , nip} C {2, 3, . . .} is a finite set, and mp = max {m : m G Tip} ■ (3.8) 

The integer mi? in (13. 8p is called the degree of the functional polynomial F. For instance, if 
^Mp = {2} or Tip = {3} the polynomial F is respectively homogeneous quadratic or cubic. 
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Every m-linear operator F^"^^ in fl3.7p is assumed to be of the form of a convolution 

F(-) (Ui, . . . , U^) (k, t)= f x^"^ (k, Ui (k') . . . (k(-) (k, fc)) (3.9) 

dk' . . . dk("^-i) 



where 



k^'™) (^k, fc) = k - k' - . . . - k(™-i), = (k', . . . , k^"^)) 



(3.10) 



indicating that the nonhnear operator F^"^^ (Ui, . . . , Vm) is translation invariant (it may be 
local or non-local). The quantities x*^™^ in (13.91) are called susceptibilities. For numerous 
examples of nonlinearities of the form similar to (13. 7p . (13.91) see [2]-[8] and references therein. 
In what follows the nonlinear term F in (12. ip will satisfy the following conditions. 



Condition 3.5 (nonlinearity) The nonlinearity F |^Uj is assumed to be of the form ([57 

/ IJ.gj) . The susceptibility x^"^^ (k, k', . . . , k*^™'-') is infinitely differentiable for all k and k.^^^ 
which are not band-crossing points, and is bounded, namely for some constant 

||x(™)|| = (27r)~('"~^)'^ sup Ix^"^ (k,k',...,k('"))| < C^, m G mt^, (3.11) 

k,k',...,k('")eK'*\o-bc 



where the norm 
defined by 



^(m) ^^k^ of m-linear tensor x*^™) : (C^'^)'" (C^'^)'" for fixed k,k 



IS 





= sup 


X(™) (k,fc) (xi,...,xj 


, where x 




|x,|<l 







Since x^™'' (k, k', . . . , k*^"*^) are smooth if k ^ (Jbc the following relation holds: 



max 

|k±k,(|<7ro, 1=1,.. .,N 



Vx[}^(k,k',...,k(™)) 



<c' 



(3.13) 



if k^i ^ (Jbc! 7i"o is defined by (13. 3p . The case when x*^"''' ^k, k^ depend on small q or, more 
generally, on q'^, q > 0, can be treated similarly, see 



3.3 Resonance invariant n/c-spectrum 

In this section, being given the dispersion relations ujn(k) > 0, n G {1, . . . , J}, we con- 
sider resonance properties of n/c-spectra S and the corresponding /c-spectra Ks as defined in 
Definition 12. 7[ i.e. 

S = {{ni,Ki), / = l,...,iV}cS = {l,...,J}xM^ ={k,,, l = l,...,\Ks\}. (3.14) 
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We precede the formal description of the resonance invariance (see Definition 13.81) with the 
following guiding physical picture. Initially at r = the wave is a multi-wavepacket composed 
of modes from a small vicinity of the n/c-spectrum S. As the wave evolves according to (12. ip 
the polynomial nonlinearity inevitably involves a larger set of modes [S]^^^ ^ 5, but not 
all modes in [S]^^^ are "equal" in developing significant amplitudes. The qualitative picture 
is that whenever certain interaction phase function (see fl4.23p below) is not zero, the fast 
time oscillations weaken effective nonlinear mode interaction and the energy transfer from 
the original modes in S to relevant modes from [5']^^^, keeping their magnitudes vanishingly 
small as /?, ^ — >■ 0. There is a smaller set of modes [SYout which can interact with modes from 
S rather effectively and develop significant amplitudes. Now, 

if [S^o^^ ^ S then 5* is called resonance invariant. (3.15) 

In simpler situations the resonance invariance conditions turns into the well-known in non- 
linear optics phase and frequency matching conditions. For instance, if 5* contains (no,k,,;g) 
and the dispersion relations allow for the second harmonic generation in another band ni so 
that 2ujnQ (k*io) = uJn^ (2k,,;(,), then for S to be resonance invariant it must contain (ni, 2k,,;(,) 
too. 

Let us turn now to the rigorous constructions. First we introduce necessary notations. Let 
m > 2 be an integer, / = (/i, .., Im), Ij € {1, ... , A^} be an integer vector from {1, . . . , A^}*" and 

C = {C^^\ , ..,C^"'^), C^^^ e {+1, -1} be a binary vector from {+1, -1}™. Note that a pair 

^ naturally labels a sample string of the length m composed of elements {^'^^\ni., k*;^ j 
from the set {+1,-1} x S. Let us introduce the sets 

A = {(C,0:^e{l,...,Ar}, Cg{+1,-1}}, (3.16) 

A'" = |a = (Ai,...,A^), G a, j = l,...,m}. 

There is a natural one-to-one correspondence between A™ and {— l,l}™'x{l,..., A^}™' and 
we will write, exploiting this correspondence 

A = ((c',/i),...,(c^"^\^^)) = (c,/), ^G{-i,ir, re{i,...,Arr for Xga- (3.17) 

Let us introduce the following linear combination 

>cm (a) = x„ (c, I) = C^'^Ki, With C^^') e {+1, -1} , (3.18) 

and let [5"]^^^^ be the set of all its values as k,,,;^ G Ks, A G A"*, namely 

|S|^.«. = U„,„,Ufc.„{-'»(^1}- (3-19) 
We call [5']^^^^^ output k-spectrum of Ks and assume that 

[5]/.,outn^bc = 0. (3.20) 
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We also define the output nk- spectrum of S by 

[^]out = {(n, k) G {1, . . . , J} X R'^ : n G {1, . . . , J} , k G [S]^^,,,} . (3.21) 
We introduce tlie following functions 

fii.m (J^ (k*^ = X^^.^-^ (^^'^^h i^*h) ' ^* = • • • ' , where k,/^, G i^s, (3.22) 

Q (^C, n, (k^^, £ j = -C^^n (k„) + Qi^rn (^A^ j , (3.23) 

where ( = ±1, m E DJlp a.s in (13.71) . We introduce these functions to apply later to phase 
functions Km . 

Now we introduce the resonance equation 

n (C, n, a) (Cx„ (a) , k) = 0, Tg {1, . . . , AT}'- , C G {-1, 1}'" , (3.24) 
denoting by P (S) the set of its solutions ^m, (, n, A j . Such a solution is called S -internal if 
(|n, (^Xj j G S*, that is n = nig, (x^ (^X^ = Ki^, Iq e {1, . . . , N} , (3.25) 

and we denote the corresponding Iq = I ^Xj . We also denote by Pint (S) C P (S) the set of 

all 5'-internal solutions to (13.241) . 

Now we consider the simplest solutions to (13.241) which play an important role. Keeping 
in mind that the string I can contain several copies of a single value I, we can recast the sum 
in (I3.22P as follows: 

n,,^ (X) = n,,^ (C, I) = , where 6, = | ^^'^'-MO ^ (|) ^ ^ , 

(3.26) 

(/) = {j G {1, . . . ,m} : = /, } , r= (/i, . . . ,U , 1 < / < AT. 



Definition 3.6 (universal solutions) We call a solution ^m, (^,72, Xj G P (S) of ([^ 



universal it has the following properties: (i) only a single coefficient out of all 6i in 



is nonzero, namely for some Jq we have 6iq = ±1 and Si = for I ^ Iq; (ii) n = nj^ and 

We denote the set of universal solutions to (13.241) by Puniv ("S"). A justification for calling 
such solution universal comes from the fact that if a solution is a universal solution solution 
for one k^, it is a solution for any other A;* G M'^. Note that a universal solution is a S'-internal 
solution with / ^Xj = Jq implying 

Puniv(5) CPi„t(5). (3.27) 
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Indeed, observe that for 5i as in fl3.26l 

implying >Cm ^Xj = ^/gk*/,, and C>Crn = 5/(,k*/(, = k,,,/,,. Then equation (13.241) is obviously 
satisfied and (^n,(xm (^jj = (^/o5k^,/p) G S. 

Example 3.7 (universal solutions) Suppose there is just a single band, i.e. J = 1, a 
symmetric dispersion relation ui (— k) = ui (k), a cubic nonlinearity F with dJlp = {3}. 
We take the nk-spectrum S = {(1, k*) , (1, — k^,)}, that is N = 2 and k^i = k*,k^2 = 

— k,f. This example is typical for two counterpropagating waves. Then ^1^3 (jc^^ = 

Zli=i C^''^^^/, (k*«,) = {Si + 62)uJi{K) and x,rn(\^ = J^'^=iC''^^Ki, = 61K1 + 62K2 = 
{61 — 62) k* where we use notation \3.2b]) . The universal solution set has the form Punw {S) = 
I ^3, C, 1, Xj : X G A^, ( = ±1 where A_|_ consists of vectors (Ai, A2, A3) of the form ((+, 1) , (— , 1) , (+, 1)) 

((+, 1) , (-, 1) , (+, 2)) , ((+, 2) , (-, 2) , (+, 1)), ((+, 2) , (-, 2) , (+, 2)), and vectors obtained 
from the listed ones by permutations of coordinates Ai, A2, A3. The solutions from Pint (S) have 
to satisfy \6i — ^2! = 1 o,nd \6i + 52I = 1 which is possible only if S162 = 0. Since ( = 61 + 62 
we have (x^ ^Xj = (5^ - k^, and (x^ ^Xj = k*i if \6i\ = 1 or (x^ (^Xj = k*2 if\52\ = 1- 

Hence Pint{S) = Punw{S) in this case. Note that if we set Si = {(l,k*)}, 5*2 = {(1,— k,,,)} 
then S = Si U S2 but Pint{S) is larger than Pint (Si) U Pint{S2). This can be interpreted as 
follows. When only modes from Si are excited, the modes from S2 remain non-excited. But 
when the both Si and S2 are excited, there is a resonance effect of Si onto S2, represented, 

for example, by X = ((+, 1) , (— , 1) , (+, 2)), which involves the mode (xm (^Xj = k*2. 

Now we are ready to define resonance invariant spectra. First, we introduce a subset 
[•5] out of [-^lout by the formula 

[S]out = {(^. k„) G [Sl^, : k„ = C^°^x„ (X) , m G DJlp, where (3.29) 
m, (,n,X] is a solution of (13.241) 



calling it resonant output spectrum of S, and then we define 

resonance selection operation 71 (S) = S U [SY^^^ . (3.30) 

Definition 3.8 (resonance invariant nA;-spectrum) The nk-spectrum S is called reso- 
nance invariant iflZ{S) = S or, equivalently, [5*]^^^ C S. The nk-spectrum S is called 
universally resonance invariant if7l{S) = S and Punw (S) = Pint (S). 
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Obviously, n/c-spectrum S is resonance invariant if and only if all solutions of fl3.24p are 
internal, that is Pint ("S") = P (S). 

It is worth noticing that even when a n/c-spectrum is not resonance invariant often it can 
be easily extended to a resonance invariant one. Namely, if TZ^ {S) fl Ubc = for all j then 
the set 

7^°° (S) = I r W (5) C S = {1, . . . , J} X 

is resonance invariant. In addition to that, TZ°° {S) is always at most countable. Usually it is 
finite i.e. 7l°° (S) = TV {S) for a finite p, see examples below; also 7?,°° (S) = S for generic 
Ks. 

Example 3.9 (resonance invariant n/c-spectra for quadratic nonlinearity) Suppose 
there is a single hand, i.e. J = 1, with a symmetric dispersion relation, and a quadratic non- 
linearity F, that is dJlp = {2}. Let us assume that k* ^ 0, k*,2k*,0 are not band-crossing 
points and look at two examples. First, suppose that 2ciJi (k^,) 7^ cui (2k*) (no second har- 
monic generation) and Ui (0) 7^ 0. Let us set the nk-spectrum to be the set Si = {(l,k^.)}, 
then Si is resonance invariant. Indeed, Ks^ = {k*}, [Si]j^ g^^f = {0, 2k=i,, — 2k*}, [SiJ^^j = 
{(1, 0) , (1, 2k*) , (1, —2k*)} and an elementary examination shows that [SiY^^^ = G Si im- 
plying 71 (Si) = Si. For the second example let us assume ui (0) 7^ and 2u!i (k*) = ui (2k*), 
that is the second harmonic generation is present. Here [SiY^^^. = {(1,2k*)} and 7l{Si) = 
{(1, k*) , (1, 2k*)} implying TZ (Si) 7^ Si and, hence. Si is not resonance invariant. Suppose 
now that 4k*, 3k* ^ a^c and cui (0) 7^ 0, cui (4k*) 7^ 2iUi (2k*), cui (3k*) 7^ coi (k*) + cui (2k*) 
and let us set S2 = {(1, k*) , (1, 2k*)}. An elementary examination shows that S2 is resonance 
invariant. Note that S2 can be obtained by iterating the resonance selection operator, namely 
5*2 = 71(71 (Si)). Note also that Puniv{S2) 7^ Pint{S2). Notice that ui (0) = is a special 
case since k = is a band- crossing point, and it requires a special treatment. 

Example 3.10 (resonance invariant nfc-spectra for cubic nonlinearity) Let us con- 
sider one-band case with symmetric dispersion relation and a cubic nonlinearity that is 
VJlp = {3}. First we take Si = {(1, k*)}, we assume that k*, 3k* are not band-crossing points, 

implying = {k*, -k*, 3k*, -3k*}. We have ^1,3 (j^ (k^^ = Z]i=i C^^^'^i (k*) = 

^icui (k*) and ^Xj = 5ik* where we use notation \3.20) . 5i takes values 1,-1,3,-3. 
// ?>uji (k*) 7^ uji (3k*) then 1^3. 24\ ) has a solution only if \6i\ = 1 and 61 = (, hence 



= k* and every solution is internal. Hence, [SiY^^^ = and 7Z{Si) = Si. Now 

consider the case associated with the third harmonic generation, namely 3uJi (k*) = ui (3k*) 
and assume that ui (3k*) + 2u!i (k*) 7^ uji (5k*), 3u!i (3k*) 7^ ui (9k*), 2uJi (3k*) + ui (k*) 7^ 
(7k*), 2c(Ji (3k*) — cui (k*) 7^ C(Ji(5k*). An elementary examination shows that the set 
S4 = {(1, 3k*) , (1, k*) , (1, — k*) (1, — 3k*)} satisfies 7Z{S 4) = S^. Consequently, a multi- 
wavepacket having 5*4 as its resonance invariant nk-spectrum involves the third harmonic 
generation and, according to Theorem \2.8[ it is preserved under nonlinear evolution. The 
above examples indicate that in simple cases the conditions on k* which can make S non- 
invariant with respect to 7Z have a form of several algebraic equations. Hence, for almost 
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all k,, such spectra S are resonance invariant. The examples also show that if we fix S and 
dispersion relations then we can include S in larger spectrum S' = TZ^ {S) using repeated 
application of the operation TZ to S, and often the resulting extended nk-spectrum S' is res- 
onance invariant. We show in the following section that nk-spectrum S with generic Ks is 
universally resonance invariant. 

Note that the concept of resonance invariant nfc-spectrum gives a mathematical descrip- 
tion of such fundamental concepts of nonlinear optics as phase matching, frequency match- 
ing, four wave interaction in cubic media and three wave interaction in quadratic media. If 
a multi-wavepacket has a resonance invariant spectrum, all these phenomena may take place 
in the internal dynamics of the multi-wavepacket, but do not lead to resonant interactions 
with continuum of all remaining modes. 



3.4 Genericity of the n/c-spectrum invariance condition 

In simpler situations, when the number of bands J and wavepackets N are not too large, the 
resonance invariance of nk- spectrum can be easily verified as above in Examples 13.91 13.1UI 
but what one can say if J or are large, or if the dispersion relations are not explicitly given? 
We show below that in properly defined non-degenerate cases a small variation of Ks makes 
S universally resonance invariant, i.e. the resonance invariance is a generic phenomenon.. 
Assume that the dispersion relations u;„ (k) > 0, G {1, . . . , J} are given. Observe then 

that Vim [c, n, = Qm ^C? (k*i, . . . , k^ii^gi) defined by (13.231) is a continuous function 
of k^.^ ^ (Tbc for every m, (, n, A. 

Definition 3.11 (cj-degenerate dispersion relations) We call dispersion relations ujn 0^) , 
n = 1,. . . , J , u-degenerate if there exists such a point k^, G M'^ \ abc that for all k in a neigh- 
borhood o/k* at least one of the following four conditions holds: (i) the relations are linearly 
dependent, namely X]n=o^"-'^" ^ ~ '^O; where all Cn are integers, one of which is nonzero, 
and the Cq is a constant; (ii) at least one of ujnik) is a linear function; (Hi) at least one of 
uOn (k) satisfies equation CuJn (k) = uon (C'k) with some n and integer C ^ ±1; (iv) at least 
one of Un (k) satisfies equation ojn (k) = ujn' (— k) where n' ^ n. 

Note that fulfillment of any of the four conditions in Definition 13.111 makes impossible 
turning some non resonance invariant sets into resonance invariant ones by a variation of k,,;. 
For instance, if 9Jti? = {2} as in Example 13.91 and 2uji (k) = uji (2k) for all k in an open 
set G then the set {(l,k*)} with k,. G G cannot be made resonance invariant by a small 
variation of k*. Below we formulate two theorems which show that if dispersion relations 
are not cj-degenerate, then a small variation of k^,; turns non resonance invariant sets into 
resonance invariant; the proofs of the theorems are given in [8J 

Theorem 3.12 // (^C, ^o, a) (^k'.^, . . . , k;,^^,) = on a cylinder G m (M^ \ (Xbc) 

which is a product of small halls Gi C (M'^ \ dbc) ^hen either ^m, C,no,Aj G Puniv{S) or 
dispersive relations u;„ (k) are u -degenerate as in Definition \3.11\ 
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Theorem 3.13 (genericity of resonance invariance) Assume that dispersive relations 
ujn (k) are continuous and not u-degenerate as in Definition \3.11\ Let /Ch™ be a set of 
points (k*i, . . . ,]<i^\Ks\) such that there exists a universally resonance invariant nk-spectrum 
S for which its k-spectrum Ks = {k,,i, . . . , k^,|7i-g| } . Then K-rinv is open and everywhere dense 

set m (M'^Vcibc)'^'''. 



4 Integrated evolution equation 

Using the variation of constants formula we recast the modal evolution equation (12. ip into 
the following equivalent integral form 

e^^^^^^^F (^U) (k, r) dr' + e'^^^'^^h (k) , r > 0. (4.1) 

Then we factor U (k, r) into the slow variable u (k, r) and the fast oscillatory term as in 
fl2:T2D . namely 

U (k, r) = e-^^('')u (k, r) , U„,c (k, r) = u„,^ (k, r) e'^^^-^, (4.2) 

where u„^^ (k, r) are the modal components of u (k, r) as in (13.41) . Notice that u„^^ (k, r) in 
(14. 2 p may depend on g and (14.21) is just a change of variables and not an assumption. 

Remark 4.1 Note that ^/Un,c (k, r) is a wavepacket, it is localized near its principal wavevec- 
tork.^. The expansion of (un (k) near the principal wavevector (k^, (we take ( = 1 for brevity) 
takes the form 

(k) = u} (k,) + Vfccu, (k,) (k - k,) + ^V^cu (k,) (k - k,)' + . . . 

To discuss the impact of the change of variables ( [^.^ we make the change of variables k— k,,, = 
^ . The change of variables Iji4-S\ ) 



U„,+ (k,r)=u„,+ (k,r)e-7-"« (4^3) 



u r(k -7-)e->"(''*)e~'?^"^"(''*)(''"^*)e"'?'>^^'^^"^''*^^''~''*^ + 



= u„,+ (k, + t r) e-^C-(k*)e-T^^-('^-)^e-T««) 
R {$) = (k) - CU„ (k,) - VkUJn {K) (k - k,) = U/lun {K) {^.f + • • • (4.4) 

has the first factor e'^^""'^^*'' responsible for fast time oscillations of U„^^ (k, r) andU„^^ (r, r) 

The second factor e~~^''^"^^*^^ is responsible for the spatial shifts of the inverse Fourier trans- 
form by ^VkUJn (k*) , since the shifts are time dependent they cause the rectilinear movement 
of the wavepacket U„,^ (r, r) with the group velocity ^VfcCUn (k*), the third factor is responsible 
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for dispersive effects. Hence the change of variables Iji4-^ effectively introduces the moving 
coordinate frame for Un,c (k, t) for every k and in this coordinate frame u„ ,^ (k, r) has zero 
group velocity and does not have high-frequency time oscillations. The following proposition 
shows that if u„_^ (k, r) is a wavepacket with a constant position, Un,+ (k, r) is a particle 
wavepacket in the sense of Definition with position which moves with a constant velocity. 

Proposition 4.2 Let ui (k, r) be for every r G [0, r*] a particle wavepacket in the sense of 
Definition \2.2\ with nk-pair {n, k^,) regularity s and position e which does not depend on 
T, assume also that and constants Ci in Ii2. 33\) and C, C in ^2.21 ) and Ii2.30\) do not depend 
on T. Let iJi (k, r) be defined in terms of ui (k, r) by Iji4-^ - Assume that 11^2. 4S ) holds. Then 
XJi (k, r) for every t G [0, r^,] a particle wavepacket in the sense of Definition \2.2\ with nk-pair 
(n, k*) regularity s and with r-dependent position r* + ^VfcCUn (k*) G M'^. 

Proof. The wavepacket (k, r) involves two components u„_^ (k, t), ( = ±1 for which 
(12:291) holds 

Un,c (k, r) = ^ {P'~y2, CK; k) n„,^ (k) u„,^ (k, r) , (4.5) 

By « 

U„,c (k,r) = u„,c (k,r)e-^^""('^). 

According to Definition (12.ip multiplication by a scalar bounded continuous function e""'''^"'-'''' 
may only change a constant C in fl2.30p . therefore it transforms wavepackets into wavepack- 
ets. To check that U; (k, r) is a particle-like wavepacket we consider (12.331) with {f3, r^,; k) 

replaced by u„^^ (k, r) e""^'^"*-*'-' and replaced by r^, + ^V^c^n (k*). We consider for brevity 
Ur, (k, r) 



r e e 



dk 



u„^^ (k, r) with C = l,the case C = — 1 is similar 
/ Vk (e'(-+^^^-"(''*))'^u„ (k, 

[ Vk ('e<'-*+f ^'=^"(i'*))'^u, (k, r) e^T^-e^^e^""^''*)) 

= / Vk(e^^*'^u„(k,r)e-^^(''-''*)) dk < h + i 
where where R (^) is defined by (14. 4p . 

h = [ e-T«(i'-i^*)Vk(e''-*''un(k,r)) dk, 

(e*-'^u„(k,r))Vke"'^^(''^''*) dk 



dk 



The integral Ji is bounded uniformly in v^hj C[P ^ ^ since u^^ij (k, r) satisfies (I2.33p . Note 
that 



h = [ (e-'^u„(k,r))Vke-^^(^-'^-) 
/" |u„(k,r)|-|Vki?(k-k,)|dk 



dk < 
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Note that according to fl4.5p and fl2.25p u„^^ (k, r) 7^ only if |k — k*| < 2/?^ and for such 
k — k* we have Taylor remainder estimate 

|Vki?(k-k,)| <C/3^-\ 



Therefore h < C [3'-' / g and 



Using (12.481) we conclude that this inequality implies (I2.33P for U; (k, r), therefore it is a 
particle-like wavepacket. ■ 

From (14. ip and (14.20 we obtain the following integrated evolution equation for u = u (k, r) , 

r > 0, 

u (k, r) = ^ (u) (k, r) + h (k) , ^ (u) = Y. •^^'"^ (^'"^ (k> ^)) , (4-6) 



^(") (u-) (k,r) = I\'t^^'^^F^ (^(e^^^-^u)™^ (k,r') dr', 



(4.7) 



where Fm are defined by (13. 7p and (13. 9p in terms of the susceptibilities x^™'\ cind JF^™-) are 
bounded as in the following lemma. 

Recall that spaces L^'*^ are defined by the formula (I2.17p . Below we formulate basic 
properties of the spaces. Recall Young's inequality 



|U* V||,^i < MIli ||v||li • 



(4i 



This inequality implies boundedness of convolution in L^'"", namely the following Lemma 
holds. 



Lemma 4.3 Let Hi,H2 G L '° be two scalar functions, a > 0. Let 



Then 



i^3(k) 



< 



Proof. We have 



sup 
k',k" 

Obviously 



Hi (k - k') H2 (k') dk'. 



^i(k) ^ i^i(k) 



< 



(4.9) 



(1 + 


k' + k"l)" 


(l + lk'l)"(l + 


k"l)" 



(l + lkl)"^ H,{k) 
(l + |k-k'|)" i^i(k-k') (1 + lk'l)'^ H2{k') 



dk'. 



1 + 


k' + k" 


1 / (1 + 


k' 


+ 


k"|) 


(1 + 


k'|)(l + | 


k"|) - (1 + 


k'l) 


(H 




k"|) 



< 1. 
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Applying Young's inequality fl4.8p we obtain 



l + |k|r ^3(k) dk< / (l + |k'|) ^i(k') dk' / (l + |k"|) ^2(k" 



dk" 



Using fl2:T8D we obtaindiS]). ■ 

Using Lemma 14.31 we derive boundedness of integral operators JF^"*) . 

Lemma 4.4 (boundness of multilinear operators) Operator T^"^^ defined hy 113. S^) . ( [^. 7| j 

is bounded from Ea = C ([0, r*] , L^'") into ([0, r*] , L^'""), a > and 



\X 



(m) I 



Proof. Notice that since L (k) is Hermitian, exp |— iL (k) ^| 
ity (US]) together with ([33]), (gT]) we obtain 

(™) (k, k 



(4.10) 
(4.11) 

1. Using the inequal- 



X 



(ui...u„,)(-,r)||^,,„ <sup 
f f [ |(l + |k'|)"ui(kO|...|(l + |k('")|)"u„(^k(") (^k,A^))|dk'...dk("^-i)dridk< 

•J M'^ •-'0 ©jTi 

||ui (ri)||^i,„ . . . ||u„ (ri)||^i,„ dn < r, \\x^'^'>\ 



R„ • • • u 



mWEa ■ 



proving (14.101) . A similar estimate produces we prove (14. lip . ■ 

The equation (14.60 can be recast as the following abstract equation in a Banach space 

u = (u) + h, u, h G Ea, (4.12) 

and it readily follows from Lemma 14.41 that JF (u) has the following properties. 

Lemma 4.5 The operator {u) defined hy {4.0^ - ^7?^ satisfies the Lipschitz condition 

||.F(ui) - ^(U2)IU, < T.Cj. ||ui - U2|U^ (4.13) 
where Cp < C^m?p {4R)^^~^ if HuiH^.^ , ||u2||^^ < 2R, with as in 1^3. 

We also will use the following form of the contraction principle. 
Lemma 4.6 (contraction principle) Consider equation 

x = J^^(x) + h, x,hG5, (4.14) 
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where B is a Banach space, T is an operator in B. Suppose that for some constants Rq > 
and < q < 1 we have 

||h|| < Ro, ||^(x)|| < i?o if ||x|| < 2Ro, (4.15) 
||^(xi)-^(x2)|| < g||xi-X2|| ||xi||,||x2|| <2i?o. (4.16) 

Then there exists a unique solution x to the equation Iji4-H\ ) such that ||x|| < 2Rq. Let 
||hi|| , ||h2|| < Rq then the two corresponding solutions xi,X2 satisfy 



^1 II 5 



|X2|| < 2Ro, ||xi - X2II < (1 - qy^ ||hi - h2|| . (4.17) 



Let Xi,X2 be the two solutions of correspondingly two equations of the form Iji4-H\ ) with T\, 
hi and T2, ^i- Assume that that J^i (u) satisfies ( [^.i5| ), ( [^.i6] ) with a Lipschitz constant 
q < 1 and that \\J-'i (x) — JF2 (x)|| < 6 for ||x|| < 2i?o- Then 

IIX1-X2II < (l-g)-'((5+ ||hi-h2||). (4.18) 

Lemma 14.51 and the contraction principle as in Lemma 14.61 imply the foUowing existence 
and uniqueness theorem. 

Theorem 4.7 Let h < -R, let < 1/Cp where Cp is a constant from Lemma \4.5\ 
Then equation lji4-6\ ) has a solution u & Ea = C ([0,r*] , //^'") which satisfies HuH^,^ < 2R, 



and such a solution is unique. Hence the solution operator u = Q (hj is defined on the ball 
< R. 

Ea 

The foUowing existence and uniqueness theorem foUows from Theorem 14. 7[ 

Theorem 4.8 Let a > 0, ^2J\) satisfy (EUP and h G L^^"" {R'^) , h < R. Then 

there exists a unique solution u to the modal evolution equation Ii2. 1\) in the functional space 
([0, r*] , L-^'°), ||u||^^ + ||9ru||^^ < _Ri (_R). The number depends on R and C^. 

Using the inequality (12.201) and applying the inverse Fourier transform we readily obtain 
the existence of an F— solution of (11.11) in ([0,t*] ,L°° (M"^)) from the existence of the 
solution of equation (12.11) in ([0,r*] ,L^). The existence of F-solutions with [a] bounded 
spatial derivatives {[a] being an integer part of a) follows from solvability in (7^ ([0, r*] , L^'"). 

Let us recast now the system (I4.6I) - (I4.7I) into modal components using the projections 
n„^(^ (k) as in (12. 9p . The first step to introduce modal susceptibilities X*'"!-''? having one- 
dimensional range in C^"' and vanishing if one of its arguments uj belongs to a (2 J — 1)- 
dimensional linear subspace in C^'' (j-th null-space of X^"^^^ ) &s follows. 
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Definition 4.9 (elementary susceptibilities) Let 

e = fn, C) G {1, . . . , Jr X {-1, ir = S'-, (n, C) G S 



(4.19) 

and x*-"^^ ^k, [ui (k') , . . . , u-m (k^"*))] 6e m-linear symmetric tensor (susceptibility) as in 
lis. 9\) . We introduce elementary susceptibilities x^"^^^ ^k, A; j : (C^"')™ C^"' as m-linear 
tensors defined for almost all k and k = (k.', . . . , k*^"^)) by the following formula 

X^^l^ (k, k) [ui (k') , . . . , u„ (k(-))] = x^^l^. (k, k) [ui (k') , . . . , u„ (k(™))] = (4.20) 

n„,c (k) x^'"^ (k, fc) [(n„,,^. (k') ui (kO , . . . , n„^_^(™) (k(-) (k, k)) u„ (k^-)))" . 

Then using fl3.5p and the elementary susceptibilities (14.201) we get 

^('") (k, k) [u, (k') , . . . , (k(-))] = T.^/::a {^^ ^) ^^'^ ' • • • ' ^™ ('^^'"^)] • 



(4.21) 

Consequently the modal components of the operators JF*^"^^ in (14.71) are m-linear oscil- 

latory integral operators defined in terms of the elementary susceptibilities (I4.2ip as follows. 

Definition 4.10 (interaction phase) Using notations from Ii3.9\) we introduce for ^ = 
n, C ) G operator 



.F;-).(ui...u„) (k,r) 




k,fc)Ii 



^(m) 



x: ;? (k, A; 



ui (k', ri) , . . . , u„ (k(™) (k, Pj , n)] d("-i)'^Mri 
wzt/i the interaction phase function defined by 

^k, j = . ^ (^k, k 

= m - C'^n. (C'k') "... - C^"^^n,„ (C^^^^k^-)) , k^™) = k^™) (k, k 

where k^™) ^k, is defined by Ii3.10\) . 

Using J^^^l in fl4:22ll we recast J^^™) (u") in the system (!M1)-(|47]) as 

.F(-) [ui . . . , u^] (k, r) = • • • ^) ' 

yielding the following system for the modal components u„^^ (k, r) as in (12. 9p 
<C (k, r) = Y.^^^^ E,^,.. (k> + Kc (k) , (r^, C) e E. 



(4.22) 



(4.23) 



(4.24) 



(4.25) 
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5 Wavepacket interaction system 



The wavepacket preservation property of the nonhnear evolutionary system in any of its 
forms (fLTI) . (I27T1) . (jMD, f!4.12p . (14.251) is not easy to see directly. It turns out though that 
dynamics of wavepackets is well described by a system in a larger space i?^^ based on the 
original equation (14.61) in the space E. We call it wavepacket interaction system, which is 
useful in three ways: (i) the wavepacket preservation is quite easy to see and verify; (ii) it can 
be used to prove the wavepacket preservation for the original nonlinear problem; (iii) it can 
be used to study more subtle properties of the original problem, such as NLS approximation. 
We start with the system (14. 6 p where h (k) is a multiwavepacket with a given nfc-spectrum 
S = {ik^i,ni) , / = !,... ,N} as in (I2.39P and fc-spectrum Ks = i = 1, . . . , \Ks\} as in 
(I2.40p . Obviously, for any / {k^i,ni) = (k^j^n/) with ii < \Ks\ and indexing ii = I for 
I < \Ks\ according to ([230]). 

When constructing the wavepacket interaction system it is convenient to have relevant 
functions to be explicitly localized about the fc-spectrum Ks of the initial data. We implement 
that by making up the following cutoff functions based on (12.241) . (I2.25P 

(k) = ^ (k, ^9k„, (3'-') = ^ (k - i9k„)) , k„eKs, t = l,...,\Ks\, ^ = ± 

(5.1) 

with e as in Definition 12.11 and (3 > small enough to satisfy 

(3^^"^ < ttq, where ttq = ttq {S) < - min dist {k^,j, abc} • (5.2) 

2 k.:,i€Ks 

In what follows we use notations from (I3.16P and 

r=(/i,...,ue{i,...,ivr, ^= (^',...,#")) G{-i,ir, a = {rj)eA^, (5.3) 

n = (m, . . . , n^) G {1, . . . , Jr , C e {-1, ir , (5.4) 
^= (n eE"" ,k= (k', . . . , k^™)) e M'", where S'" as in (ICTjl . 

Based on the above we introduce now the wavepacket interaction system 

w,,, (■) = VI/ (., ^k„j n„,^ (■) ^ (E(,,^oeA ^''') + "^"""'^ ^'^ ^' ^ ^^-^^ 

W = (Wi,+ , Wi _, W7V_+, Wat _) G E^^, W;,^^ G E, 

with ^ (■,?9k*i) ,n„,^ being as in (^^, (^M>, T defined by (USD, and the norm in E'^^ 
defined based on (12.151) by the formula 

llwll^... =^^^J|wz,^||g, i? = C([0,r,],Li). (5.6) 
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We also use the following concise form of the wave interaction system fl5.5l) 

w = J'g, (w) + h^, where (5.7) 

h = (^^ - vT/. n h \i' ^lII \i> n h) G e"^^ 

The following lemma is analogous to Lemmas 14.41 14. 5[ 

Lemma 5.1 Polynomial operator J-"^ (w) is bounded in E"^^ , (0) = 0, and it satisfies 
Lipschitz condition 

II J> (Wi) - Tq, {W2)\\^2N < Ct^ II Wi - W2||^2JV , (5.8) 

where C depends only on as in h3.11\) . on the degree of and on ||wi||^2iv + ||w2||£;2jv, 
and it does not depend on f3 and g. 



Proof. We consider every operator (w) defined by (14.221) and prove its boundedness 

< 1 and 



and the Lipschitz property as in Lemma (4.41 using the inequality exp 
inequalities (12.241) . (13. lip . Note that the integration in ti yields the factor r^, and consequent 
summation with respect to n,(,^ yields (15. 8p . ■ 

Lemma 15. H the contraction principle as in Lemma 14.61 and estimate (14. lip for the time 
derivative yield the following statement. 



Theorem 5.2 Let < R. Then there exists > and Ri [R) such that equation 

( 15. 5\) has a solution w G -E^^ which satisfies 



w 



E2N 



<RiiR) (5.9) 



and such a solution is unique. 



Lemma 5.3 Every function ,j (k, r) corresponding to the solution of ( [5. 7| j from E is 
a wavepacket with nk-pair (k^^,?!;) with the degree of regularity which can he any s > 0. 

Proof. Note that according to (15.10 and (I5.7P the function 

w;,4k,r) = ^(k,^9k„,,/?i-^)n„,,^^(k,r), ||.F(r)||^, <C, < r < r, 

involves the factor ^!i .q (k) = \E' (k — where e is as in Definition 12.11 Hence, 

n„,^. w^,^ (k, r) = if n ^ or ^9 V ^, (5. 10) 

w,,^ (k, r) = ^ (k, ^Ki,,f5'-') wi,^ (k, r) , w^,^ (k, r) = if |k - > (5.11) 

Since 

^ (k,^k„,,/?i-^) ^ {k,^Kn,p'~y2) = ^ {k,^K,„p'-') , (5.12) 
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Definition 12. II for Wi^^ is satisfied with = for any s > and C" = in (12.301) . 
Now we would like to show that if h is a multiwavepacket, then the function 



w 



k, r)= > w;^(k, r)= > (k, r 



(5.13) 



constructed based on a solution of (15. 7p is an approximate solution of equation (14.121) (see 
notation (13.161) ). We will follow the lines of [S]. We introduce 



(5.14) 



Expanding m-linear operator JF'^™) ( (X^/^ '^i,'?) ) using notations (I3.16p . (13.171) we get 

(5.15) 
(5.16) 



^"">((E„*m)'>E,,.v.,^""'(*x). where 
= wa, . . . wa„, a = (Ai, . . . , A„) G A*". 



The next statement shows that (15.131) defines an approximate solution to integrated evolution 
equation (14. 6p . 

Theorem 5.4 Let h be a multi-wavepacket with resonance invariant nk-spectrum S with 
regularity degree s, w be a solution of ( [5. 7p and w (k, r) 5e defined by Ii5.13\) . Let 

D(w) = w-J^(w) -h. (5.17) 

Then there exists /Jg > such that we have the estimate 

\\B{w)\\^<Cg + CP\ tfO<g<l, P < P,. (5.18) 

Proof. Let 

(w) = (i - ^n,^n„.'?) (w) , h- = h - ^,,^n„,^h. (5.19) 

Summation of (15.50 with respect to I, ^ yields 

^ = Zl;,^ ^i„^n„„^.F (w) + ^i„^n„,,^h. 

Hence, from (15.50 and (I5.17P we obtain 

D (w) = h- - J^^- (w) . (5.20) 
Using (12.280 and (I2.30p we consequently obtain 



< C(3' if m ^ m- 



< CP' if k„, ^ k„ 
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h- < Cip'. (5.21) 

E 

Now, to show (15.181) it is sufficient to prove that 

||^"(w)||^<C2^. (5.22) 

Obviously, 

(w) = (i - E,, ^'.^n..'^) E ^^""^ (^"^) • (5-23) 

m 

Note that 

E v^.,.n„,,. = E(„ , ^ (■' ^k*) n„,.. (5.24) 

Using (13. 5p and ( I5.14p we consequently obtain 

(i - f.-.n,.,,.) = + 5:^^^ E,„,,„.,. « (5.26) 

with S defined in (I3.14p . Let us expand now (w™) using (I5.15p . According to (I5.23p 
and (I5.26P to prove (I5.22p it is sufficient to prove that for every string A G A™ the following 
inequahties hold 

ll^oonn,^;^^'") (w^)ll < Cg^ for (n,t9) G A, and (5.27) 
||^(-,^9k,)n„,^^('") (w^)ll < CsQ, if (n,k,) G (5.28) 

We will use (I5.10p and (15. lip to obtain the above estimates. According to (I4.24p 

^(-) [w,^] (k, r) = [^A. . . . w,J (k, r) . (5.29) 

Note that according to (l510ll if Aj = 

wa, = n„,^WA,, iin = rii and = (5.30) 

Let us introduce notation 

n0 =(n,„...,n,„), ^ (a) = for A = (ij') e A"". (5.31) 

Since 

nn',tfn„,^' = 0, if n 7^ n' or 7? V ^ (5-32) 

then (I5.30p imphes 

[^Ai • • • wa„] = if ^ = (^n, c) 7^ ^ (Xj , and, hence, 

^(-)[w,](k,r) = E.,,-^"!?(A)l-A.---WAj(k,r), (5.33) 
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where we use notation (13.171) . 05.311) . Note also that 



Iln'^^J^^J'l = Oi{n' ^noT^y^C, 

,{m) 



(5.34) 



and, hence, we have nonzero I\.n','d^^^^ {^\) 



By g;22D 



(m) 



",C,€(a) 



(5.35) 



(5.36) 



WA, (k',Ti),...,wj„ (kW (k,J) ,T,)] dl"-'>''idri. 
Now we use (15.111) and notice that according to the convolution identity in (13.91) 
|WA, (k', ri)| ■ . . . ■ I WA„ (k^'") (k, fc) , n) I = if |k - ^,k,,. 

Hence the integral (15.361) is nonzero only if ( k, A; ) belongs to the set 



> m(3 



(5.37) 



= |(k,fc) : |k«-^,k,,J </3i- 



i = 1, . . . , m, 



k - "^iKu 



< m[3 



l-e 



(5.38) 



We will prove now that if (n, k*j) ^ S then for small (3 one of the following alternatives 
holds: 



either ^ (-, ^Ki) nn',^^^"J^^ (w^^) = 



or (I5.35P holds and 



k,A; 



> c> for k, A; G 5 



(5.39) 
(5.40) 



Note then since 0^ ^ ( k, j is smooth then using notation (I3.18P we get 



k,A; 



^^xv, J I (^k„, K^j \<CI3' ' for (^k, kj EBp, 

d = {'&!,..., -dm) , k„ = C ^ , T^jk^/^ = ^^9, fj , 

Hence the alternative (I5.40p holds if 



(5.41) 



(5.42) 



and, consequently, it suffices to prove that either (15.390 or (15.420 holds. Combining (I5.38P 
with (k, t^k^j) = for |k — 'd'k^^i] > P'^~'^ we find that '^i^.^J-'^"^^ [w^] can be non-zero for 
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small (3 only in a small neighborhood of a point (^^^^ ^ ['^'Ixout' ^^^^ possible 
only if 

(5.43) 



k« = (xm (^^^, = Ki e Ks- 



Let us show that the equahty 







(5.44) 



is impossible for k** as in fl5.43p and n' = n as in (15.341) . keeping in mind that (n, k^,j) ^ S. It 
follows from (13.231) and (14.231) that the equation (I5.44p has the form of the resonance equation 
(I3.24p . Since n/^-spectrum S is resonance invariant, in view of Definition 13.81 the resonance 



equation (I5.44p may have a solution only if k* 



ii, n = rii, with (n;,k*jj G S. 



Since (n, k^,j) ^ S that implies (I5.44p does not have a solution and, hence, (15.420 holds when 
(n, k*,;) ^ S. Notice that (15. 9p yields the following bounds 



< Ri, II^tWaJI^; < C. 



(5.45) 



These bounds combined with Lemma [531 proven below, imply that if (I5.42p holds then (I5.28P 
holds. Now let us turn to (KTT^ . According to (EH]) and fOTp the term oo^^' (w^) 

can be non-zero only if = k*^ ^ Ks- Since n/c-spectrum S is resonance invariant 

we conclude as above that inequality (I5.42p holds in this case as well. The fact that the set 

of all (^Xj is finite, combined with inequality (15.420 . imply (I5.40p for sufficiently small p. 

Using Lemma 15.51 as above we derive (15.270 . Hence, all terms in the expansion (15.230 are 
either zero or satisfy (I5.27P or (I5.28P implying consequently (15.220 and (I5.18p . ■ 
Here is the lemma used in the above proof. 



Lemma 5.5 Let assume that 

WA, (k',ri),...,WA^ fk("^) fk,fc),ri 



and 



k,A; 



= 0/or (k,fc) eBp 
> cj^, > for ( k. A; ) ^ Bp, where Bp as in (E 



Then 



vI/(,^'k„)H„,c<";,^(w,^) 



< 



4g 



n. 



Wa 



3 We 



+ 



(m) I 



E 



(5.46) 



(5.47) 



Proof. Notice that the oscillatory factor in (14.220 equals to 



exp <( 10 ( k. A; ) — 



id) \k,k 



exp < 10 I k. A; ) — 



ri 



Wa 



'J 1 1 _E 
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Denoting 4>nci ~ ^i,^'^n',(:X^"^} 7 = xi"^'' and integrating fl4.22p by parts with respect to ti 
we obtain 

M/ (k, ^'k„) n„, c-^™)^. (w,^) (k, r) = (5.48) 



^ (k, ^'k„) ^% \' xt^ (k, k) w,, (k', r) . . . w,,„ (k(-) (k, k) , t) d(-i)'^^ 

i0fk,A;j ^ ^ V V y y 

- / VP (k, ^'k„) / . ^ (k, k) w,, (k', 0) . . . w,,„ (k(™) (k, fc) , 0) d(--i)'^^ 
i0fk,A;j ^ ^ V V y y 

' / vi/(k,/k,o^^^^l^x^^ (k,fc)a., [w,, (kO . . . (k(-) fk,k))] ~S^~'HdT,, 

Jb i0fk,A;j ^ ^ L ^ ^ 



where -B is the set of k^*^ for which (15.381) holds. The relations (13.111) and (I2.24p imply 
X^™""* ^k, j < Using then (I5.46p . the Leibnitz formula, (15.90 and (14. 8 p we obtain 

The main result of this subsection is the next theorem which, when combined with Lemma 
15. 3[ implies the wavepacket preservation, namely that the solution u„^^ (k, r) of (14.250 is a 
multi-wavepacket for all r G [0,r^]. 



Theorem 5.6 Assume that conditions of Theorem 5^ are fulfilled. Let u„^^ (k, r) for n = ni 
and wi^^ (k, r) be the solutions to respective systems ^^^'^ ( 15. 5]) . w he defined by ( 15. 

Then for sufficiently small (3q > we have 

IK,^ - n„,,^w||^ <Cg + C'P% 0< p <Po, 1 = 1,..., N. (5.49) 

Proof. Note that u„^^ = n„ ,yu where u is a solution of (14.60 and, according to Theorem 
14. 7[ llull^; < 2R. Comparing the equations (14.60 and (15.170 , which are u = (u) + h and 
w = JF (w) + h + D (w), we find that Lemma 14.61 can be applied. Then we notice that by 
Lemma [4.5l JF has the Lipschitz constant CpT^^ for such u. Taking CpT^^ < 1 as in Theorem 
10 we obtain fICTD from Km . m 

Notice that Theorem 12.91 is a direct corollary of Theorem 15.61 and Lemma 15.31 
Analogous statement is proven in |8| for parameter-dependent equations (12. ip with 

f(u) =F(u,f?). 

The following theorem shows that any multi-wavepacket solution to ( 14.60 yields a solution 
to the wavepacket interaction system (15.50 . 

Theorem 5.7 Let u (k, r) be a solution of lli4-(>\ ) o,nd assume that u (k, r) and h (k) are 
multiwavepackets with nk -spectrum S = {{ni,'k^i) , I = 1, . . . , N} and the regularity degree s. 
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Let also = be defined by / I5. Then the functions w^^(k, r) = \l/j;^^n„j^^u (k, r) 
are a solution to the system Ii5.5\) with h (k) replaced by W (k, r) satisfying 



h (k) - h' (k, r) 



<C/?^ 0<r<r, 



(5.50) 



and if wi^^ are solutions of Ii5. 5\) with original h (k) we have the inequality 

\\wl^ (k, r) - w,,^||^, < CP', < r < r,. (5.51) 
Proof. Multiplying (14.61) by \E'i;,^n„,^^ we get 

w;^^ = ^ (-, ^9k,ij Un,,^T (u) (k, r) +^ (-, ?9k„j n„,,^h (k) , w;^^, = ^ (-, ?9k„j n„,,^u. (5.52) 

Since u (k, r) is a multiwavepacket with regularity s we have 

||u(-,r)-w'(.,r)||^, <a/3^ where w' (-, r) = ^^^^ v]/ (., ^k.^ ^ (', r) . (5.53) 

Let us recast (15.521) in the form 

w;^^ = ^ (-, ?9k„J n„,^^ (w') (k, r) + ^ (-, t9k„J n„,,, [h (k) + h" (k, r)] , (5.54) 

h"(k,r) = [.F(u)-^(w')] (k,r). 

Denoting h (k) + h" (k, r) = h' (k, r) we observe that (I5.54p has the form of (15.51) with h (k) 
replaced by h'(k, r). Inequality (I5.50p follows then from (15.531) and (14.131) . Using Lemma 
si] we obtain flCTD . ■ 



6 Reduction of wavepacket interaction system to an 
averaged interaction system 

Our goal in this section is to substitute the wavepacket interaction system (15.51) with a 
simpler averaged interaction system which describes the evolution of wavepackets with the 
same accuracy but has a simpler nonlinearity, and we follow here the approach developed 
in [8]. The reduction is a generalization of the classical averaging principle to the case 
of continuous spectrum, see [8] for a discussion and further simplification of the averaged 
interaction system. In the present paper we do not need the further simplification to a 
minimal interaction system leading to a system of NLS-type equations which is done in [8j. 



6.1 Time averaged wavepacket interaction system 

Here we modify the wavepacket interaction system ( 15. 5p . substituting its nonlinearity with 
another one obtained by the time averaging, and prove that this substitution produces a 
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small error of order q. As the first step we recast f l5.5l) in a slightly different form by using 
expansions (15.151) . fl5.29p together with (15.331) and (I5.34p and writing the nonlinearity in the 
equation (15. 5p in the form 

^ (-, ^k„J Yi^,,T r) = ^ (•' (wx) , X = {l, () (6.1) 

^:],|(A) (^a) (1^' ^) = -^"kc [-^^ • • • -^'"] r) ^ , (6.2) 



n=n{l), {■n,C) = {ni,-d) 



with jF^™-* ^ as in fl4.22l) and n [l] as in fl5.31l) . and we call JF*-™-* fwr) a decorated 

monomial ^^^^ evaluated at w^. Consequently, the wavepacket interaction system (15. 5p 
can be written in an equivalent form 

^^'^= E E ^(■'^k,0-^l™]^^^(,^)(wx)+vI/(-,^k„^^ (6.3) 

The construction of the above mentioned time averaged equation reduces to discarding certain 
terms in the original system (16. 3p . First we introduce the following sets of indices related to 
the resonance equation ( 13.24p and flm defined by (IH 



a:::,^ = {a = {r,c) eA^^-.n^ [^,ni,x) = 0} , (6.4) 

and then the time-averaged nonlinearity jF^v by 
where JF^'"'' ^.^^ are defined in (16.21) . 



Remark 6.1 Note that the nonlinearity J-'^j^^ ^ (w) can be obtained from J-'^l by an averag- 
ing formula using averaging operator At acting on polynomial functions F : (C^)^ (C^)^ 
as follows 

(AtF)^.^ = i j\-'<^o^Fj^^ (e''^i*Mi,+,e-i'^i*Mi,_, . . . , e''^- V+, e"'^" V-) dt- (6.6) 

Using this averaging we define for any polynomial nonlinearity G : (C^)^ (C^)^ averaged 
polynomial 

Ga.,,- c {u) = hm (AtG)^^^ (u) . (6.7) 

// frequencies (pj in lid. 6]) are generic, Gav.j,( (u) is always a universal nonlinearity. Note 
that J-'av,n,,i? (w) defined by l{6.5\) can be obtained by the formula ^6. 7| j where At is defined 
by formula lid. 6\) with frequencies (f>j = uJn^ (k*j^) (it may be conditionally universal if the 
frequencies (pj are subjected to a condition of the form ^6. 2S\) . see the following subsection for 
details, in particular for definitions of universal and conditionally universal nonlinearities). 
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Finally, we introduce the wave interaction system with time- averaged nonlinearity as 
follows: 

vz,^ = ^ (-, ^9k„J (v) + ^ (-, ^9k,0 n„„^h, / = 1, . . . iV, ^9 = ±. (6.8) 

Similarly to (15.71) we recast this system concisely as 

V = J-av,*(v)+h^. (6.9) 

The following lemma is analogous to Lemmas 15.11 14.51 

Lemma 6.2 Operator J^^^ q, (v) is hounded for bounded \ G E'^'^ , J^^v,^ (0) = 0. Polynomial 
operator JF^v,* (v) satisfies the Lipschitz condition 

||.^av,vl' (Vl) - J^av,* (V2)||^2iv < Ct* || Vi - V2 ||£;2]V (6.10) 

where C depends only on a in l{3.11\) . on the power of J-" and on ||vi||^2jv + ||v2||£;2]v, and, 
in particular, it does not depend on f3, g. 

From Lemma ES] and the contraction principle we obtain the following Theorem similarly 
to Theorem 15. 2[ 



Theorem 6.3 Let 



< R. Then there exists > and > such that equation 

E2N 

lid. 9\) has a solution v G E^^ satisfying ||v||^2]v < -Ri, and such a solution is unique. 

The following theorem shows that the averaged interaction system introduced above pro- 
vides a good approximation for the wave interaction system. 

Theorem 6.4 Let v^^ij (k, r) be solution of li6.S\) and w; ,9 (k, r) be the solution of Ii5.5\) . 
Then for sufficiently small [3 the v^^^ (k, r) is a wavepacket satisfying l{5.10\) . liS.ll]) with w 
replaced by v. In addition to that, there exists (3^ > such that 

||v/,^ - wi,4^ <Cg, l = l,...,N- ^ = ±, forO< g<l, < p < Po- (6-11) 

Proof. Formula (15.101) . (I5.1ip for v;^ (k, r) follow from (16. 8p . We note that w is an 
approximate solution of (16.81) . namely we have an estimate for D^v (w) = w — JF^^ ,1, — 
which is similar to (15.171) . (I5.18p : 

||Dav(w)|| = ||w- J-av,* -h||^^^ < Cf?, if < < 1, /3 < /Jq- (6.12) 

The proof of (16.120 is similar to the proof of (I5.22p with minor simplifications thanks to the 
absence of terms with \E'oo. Using (16.120 we apply Lemma and obtain (16.110 . ■ 
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6.1.1 Properties of averaged nonlinearities 

In this section we discuss elementary properties of nonlinearities obtained by formula (16.51) . 
A key property of such nonlinearities Fj^ is the following homogeneity-like property: 

F,- c (e''^^*Mi,+, e-^<^i*Mi,_, . . . , e'-^- V+. e"'-^- V-) = e'^*^*^.,C • • • , «iv+, u^-) • 

(6.13) 

The values of 0j, i = 1,. . . N for which this formula holds depend on the resonance properties 
of the set 5* which enters (16.51) through the index set ^. First, let us consider the simplest 
case when 0j are arbitrary. An example of such a nonlinearity is the function 

F2,(; (mi,+, Ml _, M2,+, M2,-) = Ml,+Ml,_M2,+- 

We call a nonlinearity which is obtained by the formula (16. 5p with a universal resonance 
invariant set S a universal nonlinearity. 

Proposition 6.5 // Fj^^ is a universal nonlinearity, then Ii6.13\) holds for arbitrary set of 
values (j)^, i = 1, . . . , N . 

Proof. Note that the definition (16.51) of the averaged nonlinearity essentially is based on 
the selection of vectors A = (^{C, I') , ■ ■ ■ , (C^"'\ ^m)) e A™ ,5 as in (16.41) . which is equivalent 
to the resonance equation (I3.24p with n = ni, ( = This equation has the form 

- CuJn (k„) + 2^^^^ 61UJ1 (Ki) = 0, (6.14) 

with 

k„ = -C 2^^^^ "^/k*/, (6.15) 
where Si are the same as in (I3.26P . If A G A™ ^ and 



WX = (WAi . . . Wa„) = I W^',;^ . . . W^{™)^^ 



then, using (16. 2p and the multilinearity of JF^*") we get 
and 

™ N 

E^% = E,=/''^^ (6.16) 

where 61 are the same as in (I3.26p . If we have a universal solution of (16. lip , all coefficients at 
every ui (k^i) cancel out {ujn (M**) also equals one of ui (k^i), namely Un (M**) = (k*/o)). 
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Using notation (13.261) we see that a universal solution is determined by the system of equations 
on binary indices 



0, / ^ Jo, 5 



lo 



jer-i(/o) 



(6.17) 



Obviously, the above condition does not involve values of uJi and Hence if 6i,(^ correspond 
to a universal solution of (16.141) then we have an identity 



(6.18) 



-.AT 



which holds for any (0^, . . . , 0^) G C 

Consider now the case where the n/c-spectrum S is resonance invariant but may be not uni- 
versal resonance invariant. Definition l3.8l of resonance invariant ?7,/c-spectrum implies that the 
set P (S) of all solutions of (13.241) coincides with the set Pint (S) of internal solutions. Hence, 
all solutions of ( 16.14p . ( 16.15^ ) are internal, in particular k^,* = k^,/p, ci;„ (k*^) = Unj^ (M*io) 
with some /q. 

If we have a non- universal internal solution of (I6.14p . ui (k*/) satisfy the following linear 
equation 

where at least one of hj is non-zero. Note that if (16.191) is satisfied, we have additional 
(non-universal) solutions of (13.241) defined by 



E c 



(i) 



5u l^Io, 5^ C 

jer-i(/o) 



(i) 



(6.20) 



Now let us briefly discuss properties of equations (I6.20p . The right-hand sides of the above 
system form a vector b = (bi, . . . ,bj\f) with 6/ = 6i, I ^ Jq, and 6/o = C + ^lo- Note that 
/ = (/i, . . . Jrn) is uniquely defined by its level sets (/). For every / the number 5+/ of 
positive (^■'^ and the number 6-i of negative (^■'^ with j E l^^ (/) in (I6.20p satisfy equations 



5. 



Si, S+i 



6. 



(0 



(6.21) 



where 



/ ^ (/) = Q is the cardinality (number of elements) of / ^(0- Hence, 6j^i,6-i are 



uniquely defined hj 6i, I ^ (/) . Hence, the set of binary solutions ( of (16.201) with a given 

b and a given I = (Zi, . . . , Im) is determined by subsets of (/) with the cardinality 6+i 
elements. Hence, every solution with a given b and a given / can be obtained from one 
solution by permutations of indices j inside every level set (/). If b is given and the 

cardinalities (/) = q are given, we can obtain different I which satisfy (]6.20p by choosing 
different decomposition of {1, . . . ,m} into subsets with given cardinalities q. For given b 
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and c = (ci, . . . , c^) we obtain this way the set (may be empty for some b, c) of all solutions 
of (16.201) . Solutions with the same h and c we call equivalent. 

When for a given wavepacket there are several non-equivalent non-universal solutions, the 
number of which is denoted by Nc, we obtain from (I6.19P a system of equations with integer 
coefficients 

^^^^ hi^,ui {Ki) = 0, ^ = 1, . . . , iV, (6.22) 
and solutions to (13.241) can be found from 

^ ((J) = 6; for some z, < i < A^'^ (6.23) 

iGr-i{/) 

where to include universal solutions we set hi^ = 0. 

Hence, when a wavepacket is universally resonance invariant, we conclude that all terms 
in satisfy fl^TTD . Since (KIM holds, we get fICT]) for arbitrary (^i, . . . , 0^) G C^. If 
the wavepacket is conditionally universal with conditions (I6.23p . then using (16.161) and (I6.23P 
we conclude that (I6.18P and (I6.13P hold if (0^, . . . , 0jv) satisfy the system the equations 

J]^^^M; = 0, ^ = l,...,iVe. (6.24) 

Now we wold like to describe a special class of solutions of averaged equations. The 
evolution equation with an averaged nonlinearity has the form 

drUj,+ = —Cj (-iV) f/,,+ + F,- + (f/i,+, f/i,_, . . . , Un,+, Un,-) , (6.25) 
Q 

drU,,. = -q (iV) U,,+ + F,, _ (t/i,+, f/i,_, . . . , Un,+, Un,-) , j = 1,...,N, 
Q 

where C (— iV) is a linear scalar differential operator with constant coefficients. The char- 
acteristic property (16.131) implies that such a system admits special solutions of the form 

Uj^(^ (r, r) = e-^-^^'/^V,- ^ (r) (6.26) 

where Vi^c^ (r) solve the time-independent nonlinear eigenvalue problem 

-i<P,V,,+ = -iC, (-iV) V;- + + + (1^1,+, . . . , Vn,-) , (6.27) 
i0^.V^,, _ = iC, (iV) V,,+ + qF^,_ (Vi,+, . . . , Vn,+, Vm,-) , j = 1,...,N. 

6.1.2 Examples of universal and conditionally universal nonlinearities 

Here we give a few examples of equations with averaged nonlinearities. When the multi- 
wavepacket is universal resonance invariant, the averaged wave interaction system involves 
NLS-type equations. 
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Example 6.6 The simplest example of l{6. 25\) for one wavepacket (N = 1) and one spatial 
dimension (d = 1) is Nonlinear Schrodinger equation 



+ 1 



(6.28) 



-a2dlUj_ + -aof/i _ + ai(9^.f/i _ + igt/i,+f/J 
Q Q 



Note that by setting y = x + air/ q we can make ai = 0. Obviously, the nonlinearity 

F^{U) = -KqUi,-^Ul^ 

satisfies ^6.13\) : 

The eigenvalue problem in this case takes the form 



-ia2dlVi^+ - iaoVi,+ + aid^Vi^+ - igqVi^^Vl 



-i0iVi _ = ia2dlVj- + iaoVi - + aid^Vi- + igqVi^+Vl_. 
If ai = and we consider real-valued Vi,4. = Vi- we obtain the equation 

(01 + ao) = -a2dlVi,+ - gqV^^ 

or, equivalently, 

^1.+ + 

gq 



gq 



0. 



c2 = > ^^^^ 



gq 



gq 



-b' < 0, 



the last equation takes the form 

with a family of classical soliton solutions 

b 



cosh (6 {x — xo) / c) 



Note that the norm of Fourier transform 



(6.29) 



(6.30) 



Cb where C is an absolute constant, 



Hence to have Vi_+ bounded in uniformly in small g according to ^6. 30\) we should take 
01 = — Oq — b'^gq with a bounded b. 
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If the universal resonance invariant multi-wavepacket involves two wavepackets (A^ = 2) 
and the nonlinearity F is cubic, that is OJli? = {3}, a semilinear system PDE with averaged 
nonlinearity has the form 

dtU2,+ = -tL2 (ZV) U2,+ + U2,+ (g2,l,+f/l,+f/l,- + Q2,2,+f/2,+f/2,-) , 
dtU2,- = lL2 i-lV) U2,- + U2,- (g2,l,-f/l,+f/l,- + Q2,2,-U2,+U2,-) , 

dtUi,+ = -tLi (^V) f/i,+ + (gi,i,+f/i,+f/i,_ + Qi,i,+U2,+U2,-) , 

9if/i,_ = (-^v) + (gi,i,-f/i,+f/i,- + gi,i,_t/2,+t/2,-) . 

Obviously, fl6.13p holds with arbitrary (j)-^, 02- 

Now let us consider quadratic nonlinearities. In particular, let us concider the one-band 
symmetric case Un (k) = ui (k) = ui (— k), i.e. J = 1, Tip = {2} and m = 2. Suppose that 
there is a multi-wavepacket involving two wavepackets with wavevectors k*i,k^.2 i.e. N = 2. 
The resonance equation fl3.24p takes now the form 

- C^i (C'k,,, + CKi,) + C'ooi {KiJ + C'uJi (Ki,) = 0, (6.31) 

where li,l2 £ {l!2}, CX'X" ^ {"1; !}• possible cases, and there are exactly four of 
them, correspond to the four well known effects in the nonlinear optics: (i) h = h, C = C" 
and (' = correspond respectively to second harmonic generation and nonlinear optical 
rectification; (ii) h h, C' = ("si-nd (' = —(" correspond respectively to sum-frequency and 
difference- frequency interactions. 

Let us suppose now that k*i,k*2 7^ and Ui (k^i) 7^ 0, Ui (k^,2) 7^ 0, where the last 
conditions exclude the optical rectification, and that k^,j 7^ and k^,j, 2k^,j, 0, ±k*i ± k^,2 
are not band-crossing points. Consider first the case when the wavepacket is universally 
resonance invariant. 

Example 6.7 Suppose there is a single band, i.e. J = 1, with a symmetric dispersion 
relation, and a quadratic nonlinearity F , that is dJlp = {2}. Let us pick two points k^,i 
and k*2 7^ ±k*i, and assume that k*i 7^ and k=Kj, 2k^,j, 0, k*i ± k*2 are not band-crossing 
points. Assume also that (i) lojx (k*j) 7^ ujx (2k*j) , Z = 1,2 , so there is no no second 
harmonic generation; (ii) cui (k*i) ± ojx (k*2) 7^ oj\ (k*i ± k*2) , (no sum/ difference- frequency 
interactions); (Hi) ujx (0) 7^ 0, ojj (k^.i) ± uji (k*2) 7^ 0. Let us set the nk-spectrum to be the 
set Si = {(1, k^,i) , (1, k^.2)}- Then Si is resonance invariant. 

In this case (16.311) does not have solutions. Hence ^ = and the averaged nonlinearity 
equals zero. 

Now let us consider the case where the wavepacket is not universal resonance invariant, 
but conditionally universal resonance invariant. In the following example a conditionally 
resonance invariant spectrum allows for second harmonic generation in the averaged system. 

Example 6.8 Suppose there is a single band, i.e. J = 1, with a symmetric dispersion 
relation, and a quadratic nonlinearity F, that is DJlp = {2}. Let us pick two points k*i and 
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k*2 such that k*2 = 2k*i, and assume that k^j 7^ and k^,j, 2k*j, 0, ±k^.i ± k^,2 are not band- 
crossing points. Assume also that (i) 2uJi (k^,i) = Ui (2k,,i) (second harmonic generation); (ii) 
uji (k*i) icuj (k*2) 7^ oJi (k*i ± k*2) , j, / = 1, 2 ('no sum- / difference- frequencies interaction) ; 
(Hi) oJi (0) 7^ 0, ujj (k*i) ± cu/ (k^,2) 7^ 0. Let lis set the nk-spectrum to be the set S = 
{(1, k*i) , (1, k*2)}- Then S is resonance invariant. The condition ^6.19\) is takes here the 
form 

2uji (k^i) - LJi (k^2) = 0, 2k^i - k*2 = 0, 



and the condition i6.24\) turns into 



luji (k*i) - uji (k*2) = 0. 
The wavepacket interaction system for such a multiwavepacket has the form 

dtU2,+ = -lL2 (ZV) f/2,+ + Q2,2,+f/l,+ ?7i,+ , 

dtU2,- = iL2 i-iV) U2,- + Q2,2,^f/l,-?7i,_, 

dtUi,^ = -tL, (^v) f/i,+ + gi,2,+f/2,+f/i,-, 

6.2 Invariance of multi-particle wavepackets 

The following Lemma shows that particle wavepackets are preserved under action of certain 
types of nonlinearities with elementary susceptibilities as in fl4.20p . In the following section 
we show in particular that universal nonlinearities are composed of such terms. 



Lemma 6.9 Let components w;. ,^ = w^- 0/ = w^^ . . . w^^ be particle-like wavepackets 
in the sense of Definition \2.Si and J^^"^^^^^^-^ {^x) Ii6.5\) . Assume that 

w,,^,c(k,/5) = \k-CKH\>P'-% C = ±, ^ = !,•••, ^- (6.32) 



Assume that vector index A G A™ ^ be such a vector which has at least one component 
\j = such that 

VuJni [Kl) = VUni. (k*;J . (6.33) 



Then for any G 



Ci 



V7 -ir,k(^) TT ,, 



< 

E 



(6.34) 



Q 



n 



where C does not depend on and small j3, g. 
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Proof. Note that 
We have by Km 



r,k = (k' + . . . + k(™)) . 



(wx) (k, r) = Vk / / exp 



icf) I k, k 



(6.35) 



^e--xi7(k,A:jw,,,,(k')...w,^^^, 



m) 



Without loss of generahty we assume that in (16.331) Ij = Im (the general case is reduced to 
this one by a renumeration of variables of integration). By the Leibnitz formula 



^,i(a: 



(w,) 



(k, r) = Ji + /2 + Is 



(6.36) 



where 



h= Vk exp <^ i0g ^ k, A; 

^X':;^ (k, fc) e— ^ w,,,, (kO . . . e— w,^,,c™) 



ir*k > X 




\Ef exp < i0g^ ^k, A; ) — - ir*k )■ x 



Vk (vl/ (k, Ku P'-') (k, fcj j J e-**^ w,,,^, (k') . . . e— • w,^_^(.) 




./ [— 7r,7r 



(m-l)d 



exp < i0g 7 k, ir^k \ x 



v&xi"^) (k, k) e— ^ w,,,, (k') . . . Vk (;e--^""'('^-^: 



Q 



k^") (k,A?) ) ) d('^-i)"'Mri 



Since ,j are bounded, we have 



(k(^')) 



< 



(k(^)) 



Using (US]) and flOTD we get 



m—l 



y -ir.k('")(k,A:)^ 

t77 



m) 



dri. 



From dSSZD, dMSD, flXT^ and the smoothness of ^ {k,Ki,P^-') we get 



(6.37) 



(6.38) 



(6.39) 
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Now let us estimate Ii. Using fl4.23p we obtain 



X 



(m) 



h 








Q 




k,k 





[— 7r,7r 



{m-l)d 



exp < 1 



k,k 



Q 



(6.40) 



-ev^uj.,, (k) + (^'"Vko;.,^ (k^'") (k, k) 



The difficulty in the estimation of the integral Ji comes from the factor y since g is small. 
Since (16.321) holds, it is sufficient to estimate Ji if. 



b-(i) _ 

According to (13.181) . since A G ,5, we have 



< (3^^' for all j. 



(6.41) 



Hence, using (I6.33P and ( 14.23^ we obtain 



0. 



(6.42) 



(6.43) 



Using (13. 2p we conclude that in a vicinity of k^, defined by (I6.4ip we have 

-^Vk^ (k) + C^"Vk^ (k(™) (k, fc))] I < 2 (m + 1) C^,2/3'"^ 
This yields the estimate 

\h\<CsP'-yg. 

Combining ([S3S]), dESSD and §M) we obtain (ICTj) . ■ 

We introduce a /5-dependent Banach space of differentiable functions of variable k by 
the formula 

l|w|Ui(r.) = P'"" II Vk (e-'^*'^w) 11^ + ||w||^ . (6.44) 
We use for 2A^- component vectors with elements Wj (k) G E"^ the following notation 

w (k) = (wi (k) , . . . , Wat (k)) , = (r*i, . . . , r^^) , (k) = (wi,+ (k) , w,; _ (k)) , (6.45) 
e--*k^ (k) = (e"--'^wi (k) , . . . , e-'^*- (k)) , 
Similarly to (15.61) we introduce the space {E^Y'' (?,) with the norm 

ll^ll(Eir(f,) = E^,^ l|w«,^lli?i(r.,) • (6-46) 
The following proposition is obtained by comparing (16.441) and (I2.33p . 
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Proposition 6.10 A multi-wavepacket w is a multi-particle one with positions r^i,. . . ,r^f^ 
if and only if 

l|w||(iJl)2^(f.) < C 

where the constant C does not depend on f3, < (3 < 1/2, and f*. 

In view of the above we will call (r*) and (E^)'^'^ (f*) particle spaces. We also use 
notations 

V^2WX = (-, . . . , VI/ (., k,,„, W,„) , 

. (W™) = ^ (-, k,;, P'-') ^('"^ .... (*2Wr) . 



Lemma 6.11 Lei w, v G (E^)"^^ (r*) ^''^^ -^^"^l ^^^^ i^x) (E^i- Assume that vector 

index X G ^ 6e snc/i a vector which has at least one component \j = (^(j, Ij) with Ij = I. 
Assume that ( (i.gj) holds and \E' (■, k*, /S^"*") is defined in ^2.25^) . Let ||w||^^i 
Then 



^ r~f II ~ ll™—! II ~ I 

< Cr* ||w||^^^2jv ||w| 



(6.47) 



where C does not depend on P, < P < 1/2, and on r*, r*;, and r^, is defined by \6.4^ . If 
II v|||.^ij2JVj.-^j < 2R the following Lipschitz inequality holds 



jrim) ^ (w) - . (v) 



< Cr* ||w 



(£;i)2^(f*) 



(6.4J 



where C does not depend on P , < P < 1/2, and on r^, r^i. 



Proof. Note that vj/gw^ and '^2^x wavepackets in the sense of Definition 12.21 To 
obtain fl6.47p we apply inequality f l6.34p and use (11.91) : for the part of E'^-norm without 
k-derivatives we use (I4.10p . Using multilinearity of ^^"^^^^ "we observe that 



(m) 



(w) - r 



m) 



n,,i?,A,<I'2 



= 5Z-^lri,A,*2 • • • '^A. - va^,^,, . . . , va„) . (6.49) 



We can apply to every term inequality (16.341) . Multiplying (16.340 by P^'^'' and using (II. 9p 
we deduce ( 16.480 . ■ 

Now we consider a system similar to (16. Sp . 



VZ,-^ = ^av,M'2,ni,i? (v) + * (", "(^krij fln,,^!!, / = 1, . . . iV, "i? = ±, 

where .Fav.^.n,,!? is defined by a formula similar to (l6.5p : 



.7^ 



(v) = E 



(m) 



The system (I6.50p can be written in the form similar to (16.90 

V = J^av,*2 (v) + h^. 



(6.50) 
(6.51) 

(6.52) 
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Theorem 6.12 (solvability in particle spaces) Let the initial data h in the averaged 
wavepacket interaction system Ii6. 52\) he a multi-particle wavepacket h {(3, k) with nk-spectrum 

S as in 3y\) . the regularity degree s and with positions r^:i, I = 1, . . . , N . Let h < 

R. Assume that S is universally resonance invariant in the sense of Definition \3. 8[ Then 
there exists t^,^, > which does not depend onv^, (3 and g such that if r^. < r** equation 
h6.52^) has a unique solution v in (E^)"^^ (r^,), such that 

\\^\\{E^f^{r,) ^ 2-R (6.53) 

where R does not depend on g, f3 and on r^, . This solution is is a multi-particle wavepacket 
with positions r^i 

Proof. Since S is universally resonance invariant every vector index A G ^ has at 
least one component Xj = (Cj, Ij) with Ij = I. Hence Lemma [6. Ill is applicable and according 
to fl6.48p the operator JFav,*2 defined by (16.511) is Lipschitz in the ball ||v||(^i)2JV(.j:^^ < 2R 
with a Lipschitz constant C't^, where C which does not depend on g, j3, and on r*. We choose 
r^.* so that C't^^ < 1/2 and use Lemma [4.61 According to this Lemma equation (I6.52p has 
a solution v which satisfies f l6.53p . This solution is is a multi-particle wavepacket according 
to Proposition I6.10[ ■ 

Theorem 6.13 (particle wavepacket approximation) Let the initial data h in the in- 
tegral equation i2.14\) with solution u(r,/?;k) he an multi-particle wavepacket h. {(3 ^\s) with 
nk-spectrum S as in Ii2. 39\) . the regularity degree s and with positions r^i I = 1, . . . , N , and 

components ofh{f3,k.) satisfy the inequality h < R. Let t^, < t^^. Assume that 

S is universally resonance invariant in the sense of Definition \3. Si We define v(r, /3;k) hy 
the formula 

N 

V (r, /5; k) = ^ T=± ^'-'^ {r,P;k), l = l,...,N, (6.54) 

1=1 

where v/^^(r,/9;k) is a solution of ( fg.gj) . Then every such Vi(k;r, /?) is a particle-like 
wavepacket with the position r^i and 

sup ||u(r,/3;k) -v(r,/3;k)||^i < Cl^? + C2/3^ (6.55) 

0<T<r, 

where the constant Ci does not depend on g, s and j3, and the constant C2 does not depend 
on g, p. 

Proof. Let v G (E^)"^^ (r^,) be the solution of equation (16.521) which exists by Theorem 
I6.12[ it is a particle-like wavepacket. Note that 
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and solution of (16.521) has the form V; ,9(r, /?;k) = \1/ k*;, [. . .] and, consequently, 
for such solutions \E'2V^ = v^, the nonlinearity •^^'"j^^ (v) coincides with and equation 
\E' (■, k,./, JF*-™^^,^. (v^) and the equation fl6.50p coincides with (16.81) . Hence, v is a 
solution of (16. 8p . Estimate (16.551) follows from estimates (16.111) and (I5.49p . ■ 

Corollary 6.14 If conditions of Theorem \2.10\ are satisfied, the statement of Theorem \2.10\ 
holds. 

Proof. Note that functions w[^(k, r) = \E'jj^^nn;,^u (k, r), = ±, in Theorem 15.71 are 
the two components of u/ (r,/5; k) in (12.451) . Hence, (I5.5ip implies that 

||u; - w,,+ - w;,_||^ < C"/3^ < /3 < /3o, (6.56) 

where w; ,? are solutions to (15. 5p . According to (16.111) . if v/ ,j (k, r) is the solution of (16. 8p 
we have 

II V - ^i,4e <Cq, l=l,...,N- ^ = ±. (6.57) 

Hence, 

\\xii-wi^+-^i^^\\^<CQ + C'P% 0</3</?o- (6.58) 

This inequality implies (I2.46p . We have proved that v; ,^ is a particle-like wavepacket as in 
Theorem 16.131 Estimate (16.581) implies that u; is equivalent to v; = v; _|_ + v; _ in the sense 
of (I2.42p of the degree Si = min (s, Sq). m 

7 Superposition principle and decoupling of the wavepacket 
interaction system 

In this section we give the proof of the superposition principle of [9] which is based on a 
study of the wavepacket interaction system (16. 8p . We will show that when we omit cross- 
terms in the averaged system wavepacket interaction system, the resulting error is estimated 
by |ln/3|, that is component wavepackets evolve essentially independently and the time 
averaged wavepacket interaction system almost decouples. 

Let JFav,n,,i9 be defined by (16.50 and let a decoupled nonlinearity JFav,ni,i9,diag be defined 

by 

^av,n„^,diag (w) = Y^rnem^ "^S' -^S.diag (w) = J^XeA^^^T' ■^"^ A«(^) ^^^^ ' '^^ 

where the set of indices K^f^^^ consists of 

a:::^^ = {^={l c) e a;^:, :/, = /, j = i, m} . (7.2) 
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Note that ^^^^^lag (17.11) depends only on w/^+ and Wi^_: 

•^S.diag (w) = ^Sag,/ (W,) , Wi = (W«,+, W^,.) . (7.3) 

The couphng between different variables v; in (16. 8p is caused by non-diagonal terms 

) • (7.4) 

Obviously, equation (16.91) can be written in the form 

(7.5) 

The system of decoupled equations has the form 

Vdiag = ^av,*,diag (Vdiag) + . (7.6) 

or, when written in components, 

Vdiag,^ = ^S,diag,« (Vdiag.O + h^,, / = 1, . . . , A^. (7.7) 

We will prove that contribution of J^av,*,coup in (17.51) is small, the proof is based on the 
following lemma. 

Lemma 7.1 (small coupling terms) Let •^^"ji'^^^ ('^a) ( f^. 5\) . let all compo- 

nents of satisfy i[6.32\} and be wavepackets in the sense of Definition \2.1\ let also 
/ li.gj) hold. Assume also that: (i) the vector index A has at least two components Aj = {Ci,h) 
and Xj = (CjJj) with li ^ Ij; (ii) both w^. and are particle wavepackets in the sense of 
Definition \2.2\ ' (Hi) either ^2. 51\) or { 2.54 ) holds. Then for small (3 and g 



n,,i9,5(A) V X/ 



<cj-\\np\. {7.1 



EN 



p 



Proof. Since k*; are not band-crossing points, according to Definition [STT] and Condition 
13.21 the inequalities (13. 2p and (I3.13P hold. According to the assumption of the theorem at 
least two w;^, are different for different j. Let us assume that Ij-^ = li, Ij^ = Im , h (the 
general case can be easily reduced to this one by a relabeling of variables). Since w;^and w^^ 
are particle wavepackets, they satisfy (12.330 with r replaced by r^^ and r^^ respectively. 
Let us rewrite the integral with respect to ri in (I4.22p as 

dA)(-x)(k>-)= (7.9) 

L L {^^«^ 7} ''"^''"''^^ 

(k, k) = x[J (k, k) w,, (k') . . . (k(™)) , (7.10) 
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where 



and then rewrite fl7.9p in the form 



where 



/ / exp^(k,i,Ti,e,r,,,r,„~)^(k,i,r,,,r,„~)d<"'-i>''&Ti. 

exp^ (k, k, Ti, r;,, r;^ j = exp |i0^_^ (^k, ^) ^ - i^^ik' - ir^^k^™)! , 
A (k, fc, r,,, r,„) = e^r^.k'^ir.^kC-.^M j^j^^ _ 



(7.11) 



(7.12) 



According to (13.101) k*^™^ ^k, = k — k' — . . . — k^™ . Hence, picking a vector p with a 
unit length we obtain the formula 



exp^ (k, A;,ri,^,r;^,r;„ 



■ Vk' exp^ ( k, A;, r 1 , ^, r;^ , r^,, 



P ■ Vk'0^_^ ( k, ) Ti - ^p ■ (r^, - 



(7.13) 



If we set 



6*0 k, /c, ^,ri 



then (17.131) can be recast as 



0' = Vk'0^,e (k,z, = Vk'cu (C'k'J - Vk(™)a; (C^"^^kl-) 
Cp = p ■ 0', gp = ^p ■ (r/, - rz„) , 
(cpTi - gp) 



P • Vk'0^_^ ( k, ) ri - p ■ (r;, - 



(7.14) 



(7.15) 



^p • Vk' exp^ f k, k, Ti, vi^.vi 
exp^ ( k, k, ri, q, r^,, r/„ ) = ry— -y ^^o ( k, /c, q, Ti 



i (CpTi - Qp) 

If (I2.5ip holds the vector 0' 7^ and to get |cp| 7^ we can take 

P = I0'r^ ■ 0'. |cp| = Po > 0. 
If (I2.54P holds we have (f)' = and we set 



1-1 



(7.16) 



(7.17) 



(7.1^ 



Let consider first the case when (I2.5ip holds. Notice that the denominator in (I7.16P vanishes 
for 

rio = ^. (7.19) 
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We split the integral with respect to Ti in fl7.11l) into a sum of two integrals, namely 



Fi= [ [ expJk,k,T,,g,ri,,ri]A(kJ,ri,,ri]~S^-'^''kdT,, (7.20) 

where cq is a large enough constant which we estimate below in (17.281) . Since wj are bounded 
in E and (I2.48P holds, we obtain similarly to (I4.10p the estimate 



< Cco(3'~^ \ln(3\ n < Ci (R) 



f7.211 



To estimate the norm of Fi we use (I7.13P and integrate by parts the integral in (I7.20p 
with respect to k'. We obtain 



/(k,ri) 



Fi= / /(k,ri)rfri, 

■>'|rio-ri|>/3i-'|ln/3| 

^exp^ (k, k,Ti,g, ri^,ri 



(7.22) 



1 [CpTi - QQp) 



P - Vi. 



9oA [k,k,ri^,ri 



According to (I7.10p and (I3.10p the expansion of the gradient Vk' in the above formula 
involves derivatives of x, Oq, e^^^^^'wi-^ and e"^'™-'''™'w;^. To estimate Oq and V^o we note that 



k, k, g,Ti 



[p ■ (/)Vi - qp) 



(p ■ (/)Vi - qp) + rip ■ Vk'0^_^ (k, kj - cf)' 
1 



(7.23) 



l + Tip- (^k, A;j - 0' /{cpTi-qp) 

Since |rio — ti| > cqIS^'" |ln/5|, from (17.191) we infer 

\cpTi - qp\> CpCoP^'^ |ln/3| . 
From (I6.32p . we see that in the integral (17.220 integrands are nonzero only if 



s * 



< TTo/S^-Mk-Ck^l <m7ro/9 



l-e 



(7.24) 



(7.25) 



where ttq < 1. Using the Taylor remainder estimate for Vk'0^ ^ at A;* we obtain the inequality 



Vk'0(._^ {k,k] - (j)' 



< 2ma,2/5'-^ 



(7.26) 
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Hence in (Ol 



Tip ■ 



Vk'0(^^^ (k, k) - <p' I (cpTi - gp) < 2mT^C^^2l (cpCo |ln/?|) . (7.27) 



(7.28) 



Suppose that [3 < 1/2 is small and Cq satisfies 



|ln/3| - ln2 - A^"'"^''"- 



Then it follows from fl7:23|) with help of ([T^HD, (O, (17^ and (17:271) that 



(k, k, g,Ti 



< 2. 



(7.29) 



Obviously, 

Vk'6'o ( k, k, g, ti 



-TlV, 



p- ( ( k, k) - (j)' 



(CpTi -Qp) 1 + Tip- 

Using (II27D, (17:281) and (Q we obtain that 

4r. 



( k, k) -4)' I (cpTi - gp) 



(7.30) 



Vk'6'o k, k, g, ti 



< 



\CpTi 



P ■ ( Vk'0^^f (k,/c) - 0' 



< 



.. (7.31) 



|cpri - qp\ 



To estimate Vk'X we use ( 13.131) . We conclude that the absolute value of the integral (I7.22p 
is not greater than 



|/(k,ri)|< 
2gT^ 



riCp 



A k, A;,ri,,r;,„ 



< 



ViCp - qp\ 
4C^,2gr.^ ||x^"Uk, )" ■ 



Vk'^ k, A;,ri^,ri^ 



(j(m-l)d^ 



(7.32) 
(7.33) 



FlCp — 

2^r. Ilx^") (k, 



kiCp — ^pl 

hat ||wj||j^ 
(I2.33p . Hence we obtain 



m 

b=2 



Vk'e'^'i'^'wi, 



|(V,.-V,(™))x('") (k,-)| 



m— 1 

n 

i=i 



w 



V7 ir; kf™) 

Vk(™)e '™ 



Note that ||wj||j^i are bounded according to (12.271) and Vk(m)e^'"''"^*'"'wi^, Vk'e''"'i'^'w;^ by 



|/(k,ri)| < + 



TiCp-qp \T^Cp-qp[ 



(7.34) 
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Obviously, 



I 1-1 , 1 r dn 1 , T,-qp/cp 

ri-9p/cp[>co/3i-^|ln/3| Cp Jco/jl-^lln/?! ^1 Cp Co/3 '|ln/3| 

<-{C+ |ln [/^^-^ |ln/?|] I) < - [C + \\np\ + |ln |ln/?||] < - [C + 2 |ln/?|] . 



Similarly, using fl2.48l) we get 



\Ti-qp/cp\>coP'-''\lnP\ 



,-2 , 1 r-"^^/'*' rfri 1 
riCp - Qpl dTi = — / — = — 



< 



1 



< 



Hence, we obtain for small (3 



T^^^ (wi...wj (k,r) 



E - ^13^+' 



|ln/3|. 



(7.35) 



Now let us consider the case when f l2.54p holds, (p' = and p is defined by f l7.18p . Turning 
to expression (17.231) we notice that 



CpTi 



-Q \rh -rij, \CpTi -qp\ ^ < 



l+e' 



and, according to (I7.26p . 



Then we estimate the denominator in (I7.23P and (I7.30p using (12.541) : 



np ■ 



( k, k 



<t>' /{cpT,-qp) <r,C^,2l3'-'/{Q\vi,-viJ)< 



If (3 is so small that (17.280 holds we again get (I7.29P and (I7.3ip . Hence, we obtain (I7.35P 
in this case as well (in fact, in this case the logarithmic factor can be omitted). Finally, we 
obtain (I7.36P from (I7.35P after summing up over all A, C- ■ 

Lemma 7.2 Let the nk-spectrum S be universally resonance invariant. Let the operators 
^av,n,,^> (w), J^av,n,,i»,diag (v) and J^s.v,ni,-d,coup be defined respectively by (EJj, (g7^ and ( f^TD - 
Let V, ||v||^jv < 2i?, be a multi-wavepacket solution of ^6. 9^ with the nk-spectrum S. Then 



for small (3 and g 



|^av,n,Acoup (v)||^iv < ^ l^^/^l • 



(7.36) 
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Proof. According to fl6.5p . fl7.1l) and fl7.4p J?-'av,n;,i9,coup involves only terms with A G 
A™'^ \ ^ni'^^^ is sufficient to prove the estimate (17.81) for indices A G A™'^ \ A™'^"^^. 

Such indices involve at least two components Aj = (Cj, /i) and \j = (Cj, ^j) with li ^ Ij since 
the n/c-spectrum is universally invariant, see (13.261) . According to Theorem 16. 131 the solution 
V is a particle-like wavepacket, therefore all components of are particle-like; (16.321) holds 
according to (15.111) . Hence, all conditions of Lemma [7. II are fulffiled and (17.361) follows from 

dUD. ■ 

Note now that every equation (17. 7p is an approximation of the equation (14.61) with single- 
wavepacket initial data h^, namely 



Ui{k,T)=J^{ui)(k,r) + hi (k). 



(7.37) 



One can apply to this equation Theorems 15.61 and 16.41 formally restricted to the case = 1 
of a single wavepacket. Based on this observation and on the above Lemma we prove the 
following theorem which implies previously formulated Theorems 12.141 and 12.151 

Theorem 7.3 Assume that the multiwavepacketh = ^ h; zs particle-like and its nk-spectrum 
is universally resonance invariant. Assume also that either Ii2.51\) or l[2. 54\ ) holds. Let u be 
solution of equation (^7^. Let be solutions of ^7.31 ). Then the superposition principle 
holds, namely 



u 



N 

1=1 



< c- 



1+e 



(7.38) 



Proof. Let vjiag,; be a solution of the decoupled system (17.71) . We compare systems 
rsj) and (17. 6p . The difference between the systems is the term J^-kv.nj.iS.coup (v) . According 
to Theorem 16.121 the solution v is a particle-like wavepacket and we can apply Lemma I7.2[ 
According to this Lemma ( I7.36P holds. Applying Lemma IT6] to the equations (I7.5P and (17. 6p 
and using (17.361) we conclude that the difference of their solutions satisfies the inequality 



|Vi — Vdiag,i|| E 



< c 



\\np\+C'l3\ 



(7.39) 



According to Theorem 16.131 inequality (16.550 holds where v is a solution of (16. 9p which can 
be rewritten in the form of (17. 5p . From (I6.55P and (I7.39P we infer 



TV 



Vdiag,/ 



1=1 



/3 



|ln/3| +Ci/3^ 



(7.40) 



Note that equation (17. 7p for Vdiag,/ coincides with the averaged equation (16.91) obtained 
for the wave interaction system derived for (17.371) . Therefore, applying Theorems 15.61 and 16.41 
to the case = 1 and h = h/ we deduce from (15.490 and (16.110 the estimate 



Finally, from (ElQD and ((7311) we infer (ESHl). ■ 



(7.41) 
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7.1 Generalizations 



In this section we show that the particle-hke wavepacket invariance can be extended to the 
case when rafc-spectra 5* are not universally resonance invariant. So suppose that n/c-spectrum 
S is resonance invariant and consider nonlinearities of the form similar to (16.51) 

(w) = ^S. = E^.. _ , (*^^) . (^-^2) 

where AJ^^ ,j C A™^ is a given subset of A™. Obviously, jF^v defined by ( 16. Sp has the form of 
(I7.42P with A^^ ,5 = A™ ,5. Let us introduce a multi-wavepacket 

w = (w„j,+, w„,,_, . . . , w„^,+ , w„^,_) (7.43) 

with the nfc-spectrum S = {{ni, 9) , I = 1, . . . , N; 9 = ±}. 

We call a subset S' G S sign-invariant if with it has {ni, 9) as an element then (n/, —9) is 
also its element. Suppose that S" C S* is sign- invariant. It is easy to see that if a set S" C S" is 
sign-invariant then it is uniquely defined by a subset of indices /' = /' (S") C / = {1, . . . , A^}, 
namely 

S' = {{ni,9):lel'{S'), 9 = ±}. 

Definition 7.4 We call index pair {ni,k^i) Group Velocity Matched (GVM) with Tres,nu-d 
if every nonzero term J^*'™^ ^^^^^ in the sum ^1.4^ has the index X such that for at least one 

component Xj = (^C^''\l'j^ of this index 

VuJn, {Kl) = VuJn,^ [Kl,) ■ (7.44) 

We call S' a GVM set with respect to the nonlinearity J-'res defined by (TJ^ if S' G S is 
sign-invariant and every (n;,k*i) G S' is GVM. 

Obviously, if S is universally resonance invariant and ^'^^^ = &s in (16. 5p then S is 
a GVM set, and in this case Ij = Jq as in Definition 13.61 U S' G S is sign-invariant we call 
a multi-wavepacket w as in (17.431) with the nfc-spectrum S = {(n;, 6^) , / = 1, . . . , A^; 6* = ±} 
partially S' -localized multi-wavepacket if for every {ni,9) G S' the wavepacket w^^ e is a 
spatially localized with a position r*/. Note that according to Definition 12.71 if S" = 5* is a 
partially S"- localized multi-wavepacket then it is a multi-particle wavepacket. 

Theorem 12.101 on the particle-like wavepacket preservation can be generalized as follows. 

Theorem 7.5 (preservation of spatially localized wavepackets) Assume that conditions 
of Theorem \2.9\ hold, in particular the initial datum h = h (/3, k) is a multi-wavepacket 
with nk-spectrum S. Assume also that S' G S is a GVM set, h = h (/?, k) is partially S'- 
localized wavepacket with positions r*^, / G I' {S'), and that holds. Then the solution 



u(r, /?) = Q ( jF(p(/?)) ,h(/5) j (r) to i2.14\) for any t G [0, r*] is a multi-wavepacket with 
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nk-spectrum S, and it is an S' -localized wavepacket with positions r^.^, I G /'(S"). Namely, 
Iji2.4(>] ) holds where ui is wavepacket with nk-pair {ni, k*;) G S' defined by Iji2.45\ ), the constants 
C,Ci,C2 do not depend on r*;. and every ui, I G I' {S'), is equivalent in the sense of the 
equivalence \2.42\} of degree si = min (s, sq) to a spatially localized wavepacket with position 

Proof. The proof of the Theorem is the same as the proof of Theorem 12. 101 since it used 
only the fact that a universally resonance invariant set is a GVM one, that allows to apply 
Lemma 16.91 One also have to use the space (E^) (?„ S') with the norm defined by the 
formula similar to (I6.46P : 

\MiE^ri.„s, = E,. + P'-"' E.,,,,, Il^'' (e-^^'^'^w,,) . (7.45) 

After replacing (E^)'^^ (r*) with {E^)^^ (r^,, S') we can literally repeat all the steps of the 
proof of Theorem 12.101 and obtain the statement of Theorem 17.51 ■ 

Below we prove that the superposition principle can hold not only for universal resonance 
invariant multiwavepackets, but for other cases allowing resonant processes such as second and 
third harmonic generations, three-wave interaction etc. Here we prove a theorem applicable 
to such situations, which is more general than Theorem I2.14[ 

Let us consider a multi-wavepacket with resonance invariant n/c-spectrum 

S={{nuKi), / = 1,...,A^} 

as in (I3.14p . and assume that is a union of spectra Sf 

S = SiU...USk, Spf}Sq = if p^q. (7.46) 

Recall that resonance interactions are defined in terms of vectors A G A™ (see (13.161) . (13.171) ). 
We call a vector A = ^(^', Zi) , . . . , {c,^^\^m^^ £ A™ cross-interacting (CI) if there exist 

at least two indices [c'"^\li^ and such that ^C^*\ j G Sp., {c'^^\l'j^ & Sp^ with 

Pi T^Pj- 

Definition 7.6 (partially GVM decomposition) We call decomposition l[7.40) partially 
GVM with respect to Tres defined by { 7-4^ if the following two conditions are satisfied: (i) 
every spectrum Sj, j = 1, . . . , K, is resonance invariant; (ii) a solution ^m, n,)Cj G P (S) 

of the resonance equation ![3.24^ with CI vector A = ^(C', /i) , . . . , {^^^\ j has at least two 

indices ^C*"*\ ^ and ^C*'''\ j ^ Sp. with pi ^ pj such that both li and Ij are CVM with 

respect to Tres dlT'd 

VkUn,. {KiJ - V^Un,. [Ki^] ^ 0. (7.47) 
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Now we use Lemma 17.11 for small coupling. Being given a partially GVM decomposition 
f l7.46p we introduce the set of coupling terms between Sp^ and S'p^.as follows: 



A m,coup 



I A = (/■; C) e : 3 z ^ J such that G S^^, I, G S^^ } , (7.48) 



We also introduce a set of interactions reducible to every Sp (block-diagonal) which is similar 
to (D: 

. m,red _ .m \ a m, coup (7AQ\ 

and the reduced operator 

^av,n,.,red (w) = Y^m^mi, -^S.red (w) , -^3,,,^ (w) = Y^XeA^'X ^-^^^ ^^^^ ' ^^'^^^ 

where A™'^*^*^ is defined by f l7.49l) . Note that if the set 5* is universal resonance invariant 

and every Sp. is a two-point set {(+, U) , (+, k)} then A™'^^'''^ = A^^'^^^. We introduce also a 
partially decoupled, reduced system similar to (17.61) 

Vred = -^av,*,red (Vred) + , (7-51) 

which can be rewritten in the decoupled form, similar to (17. 7p : 

Vred,p = ^S,red,p (Vred,p) + h ed,*,p> P = I, ■ ■ ■ K. (7.52) 

Now Vred,p may include more than one wavepacket, namely 

ini,e)eSp {ni,e)eSp 
The following theorem is a generalization of the Theorem 12.141 on the superposition. 



Theorem 7.7 (general superposition principle ) Suppose that the initial data h of \2. 14\ ) 
is a multi-wavepacket of the form 

K 

h = ^ h,ed,p, (7.54) 

p=l 

where h. is a multi-wavepacket in the sense of Definition \ 3.^ with resonance invariant nk- 
spectrum S, hred,p is a multi-wavepacket with a resonance invariant nk-spectrum Sp and 
the decomposition Iji7.46\ ) is partially GVM in the sense of Definition 7.6 with respect to 
the nonlinearity Tav defined by Ii6. 5\) . Suppose also that Iji2.48\ ) holds. Then the solution u 

= Q ^h^ to the evolution equation \2.14^ satisfies the approximate superposition principle: 

K 



g J2 h-d,p ]=Y^ (^-d,p) + D, (7.55) 



vp=l / p=l 
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with a small remainder D (r) satisfying the following estimate 



sup 

0<r<r, 



DM 



|ln/3| 



(7.56) 



where e is the same as in Definition \2.1\ and can he arbitrary small, does not depend on 
P, Q and e. 



Proof. The proof of Theorem 17.71 is similar the proof of Theorem 17.31 Averaged system 
(16. 9p can be written similarly to (17. 5p in the form 



^ — T =iv,*,rcd (v) + T a,v,*,coup (v) + . 



(7.57) 



Comparing now systems (I7.57P and (I7.5ip we find that the difference between them is the term 
^av,n,,i?,coup (v) . According to Theorem 17.51 the solution v is a spatially locahzed wavepacket 
and, hence, we can apply Lemma [7.21 getting the inequality (I7.36p . Applying Lemma [4.61 to 
the equations (17.570 and (I7.5ip and using (17.360 we conclude that the difference of their 
solutions satisfies the inequality 



(7.58) 



According to Theorem 16.131 inequality (16.550 holds where v is a solution of (17.570 and we 
infer from ( I7.58P 



K 



U — Vred.p 



<C^^J\nP\+C^P' 



(7.59) 



Similarly to (I7.37P we introduce equation for u 



rcd,p 



^red,i 



Ured.p (k, r) = J" (Urcd,p) (k, t) + hrcd.p (k) . 

Applying Theorems 15.61 and 16.41 we infer similarly to (I7.4ip the inequality 



(7.60) 



(7.61) 



Finally, from (173^]) and ([7SID we infer f[73BD . ■ 
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